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EXECUTIVE  SUMMARY 

One  of  the  factors  affecting  future  water  supply  is  the  nature  of  the  future  climate  regime.  The  effects  of  climate 
change  on  water  yield  from  the  basins  of  Alberta  will  affect  water  uses  and  water  management  in  those  basins. 
Alberta  is  embarking  on  a  regional  planning  process  and  the  effect  of  climate  change  on  watershed  responses  is 
a  key  consideration  in  that  process.  Hydrologic  models  have  commonly  been  used  to  assess  the  effects  of 
climate  change  on  watershed  hydrology.  As  part  of  a  regionally  specific  initiative  for  the  Lower  Athabasca 
Regional  Plan,  Alberta  Environment  (AENV)  has  identified  the  need  for  a  hydrologic  model  to  assist  it  in 
assessing  the  effects  of  climate  change  on  water  yield  in  the  Lower  Athabasca  Regional  Plan  Area  (LARP),  and 
in  the  Athabasca  River  Basin  and  the  Beaver  River  Basin,  in  particular.  The  intent  is  that  the  outcome  of  this 
regional  initiative  can  be  used  as  a  prototype  to  support  similar  projects  in  other  regions.  Further  to  a  Request 
for  Proposals  and  the  submission  of  proposal  by  Golder  Associates  Ltd.  (Golder),  AENV  contracted  Golder 
Associates  Ltd.  (Golder)  to  develop  criteria  for  the  selection  of  an  appropriate  hydrologic  model,  to  develop  an 
approach  to  implement  the  selected  model,  to  apply  the  selected  model  to  the  Athabasca  River  and  Beaver 
River  Basins  using  five  climate  change  scenarios,  and  to  assess  the  results  of  the  simulations. 

The  model  selection  criteria  included  considerations  of  a  number  of  factors  including  the  (1)  range  and 
mathematical  representation  of  hydrologic  processes  in  the  model  and  the  temporal  scale  of  the  simulations,  (2) 
applicability  of  calibrated  model  parameters  or  current  simulation  modules  even  under  modified  climate  regimes, 
(3)  flexibility  of  model  to  simulate  hydrologic  processes  over  a  wide  range  of  modeling  scale;  (4)  extent  of 
climate,  watershed,  vegetation,  and  soil  data  requirements  for  running  the  model,  (5)  ability  of  input  data 
structure  of  model  to  support  the  format  of  downscaled  climate  change  scenario  data;  (6)  outputs  of  model  and 
their  relevance  to  support  water  and  watershed  management  in  Alberta;  and,  (7)  currency  or  peer-acceptance  of 
model.  Ten  hydrologic  models  were  evaluated  for  the  purposes  of  this  study:  WATFLOOD,  MISBA,  MIKE-SHE, 
HSPF,  HEC-HMS,  HYDROTEL,  SSARR,  DHM-RS,  VIC  and  CHRM. 

The  three  top-ranked  models  were  WATFLOOD,  MIKE-SHE  and  HSPF.  MIKE-SHE  is  ranked  first  among  the 
three  because  of  its  physically-based  approaches  for  simulating  hydrologic  processes,  sophisticated  hydraulic 
routing  capabilities,  and  excellent  technical  support  and  documentation.  The  scores  for  WATFLOOD  and  HSPF 
were  not  too  different  from  each  other.  Both  MIKE-SHE  and  WATFLOOD  are  proprietary  softwares  with  no 
option  for  code  modification  by  the  user.  HSPF  is  the  other  model  in  the  top-tier.  The  representation  of 
hydrologic  processes  in  HSPF  is  not  as  sophisticated  as  those  in  MIKE-SHE,  however,  they  are  still  fairly 
rigorous  and  are  widely  acceptable  and  used.  When  only  selection  criteria  1  ( Criteria  relevant  to  objectives  of 
study  -  hydrologic  effects  of  climate  change)  and  3  ( Criteria  Determining  the  Degree  of  Sophistication  of  the 
Processes  Simulated  in  the  Model)  are  considered,  HSPF  and  WATFLOOD  are  scored  very  close  to  each  other, 
with  MIKE-SHE  scoring  higher  than  either.  When  only  selection  criterion  3  is  considered  on  its  own,  the  three 
models  are  ranked  as  follows:  1:  WATFLOOD,  2:  HSPF,  and  3:  MIKE-SHE.  HSPF  is  a  public  domain  model. 
Since  its  initial  development  nearly  twenty  years  ago,  the  HSPF  model  has  been  applied  throughout  North 
America  and  numerous  countries  and  in  various  climatic  regimes  around  the  world.  The  modules  for  simulating 
hydrologic  processes  within  the  model  (e.g.,  snowmelt  runoff,  infiltration,  etc.)  are  reasonably  sophisticated  for 
the  purposes  of  this  project.  Based  on  the  model  selection  criteria  and  further  evaluations  of  the  top-ranked 
models  and  full  weighting  of  considerations  for  the  purposes  of  this  particular  project,  and  discussions  with 
AENV,  the  HSPF  model  was  selected  for  implementation  on  the  Athabasca  River  and  Beaver  River  because  of 
its  overall  good  compromise  between  sophistication,  user-friendly  attributes,  and  wide  use.  MIKE-SHE  and 
WATFLOOD  were  the  next  most  appropriate  models  in  order  of  preference. 

The  HSPF  model  was  calibrated  and  validated  for  application  to  the  Athabasca  River  Basin  and  Beaver  River 
Basin.  Each  basin  was  sub-divided  on  the  basis  of  drainage  network,  locations  of  gauging  stations,  and  surficial 
geology.  Temperature  and  precipitation  data  for  the  calibration  of  the  model  were  from  climate  stations  closest 
(specific  selection  based  on  availability  of  concurrent  climate  and  flow  data)  to  the  sub-basins.  Seven  sub¬ 
basins  (tributary  to  the  Athabasca  River)  were  selected  for  calibration  of  the  model.  Six  nodes  located  on  the 
main  stem  of  the  Athabasca  River  were  selected  for  validation  of  the  calibrated  parameters.  The  calibration  on 
mean  annual  and/or  mean  open-water  flows  is  generally  good  for  most  calibration  sub-basins.  With  respect  to 
mean  monthly  flows,  the  calibration  is  generally  reasonable,  although  there  are  some  significant  differences 
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between  observed  and  simulated  values  for  some  winter  months.  The  calibrated  model  reproduced  the 
measured  discharges  at  the  validation  nodes  on  the  main  stem  of  the  Athabasca  River  well  to  reasonably  well, 
giving  confidence  in  the  use  of  the  model  for  assessing  the  hydrologic  effects  of  potential  future  climate  changes, 
it  was  concluded  that  the  calibrated  HSPF  model  has  been  validated  and  is  appropriate  for  assessing  the  effects 
of  climate  change  on  the  yield  in  the  Lower  Athabasca  Regional  Plan  Area. 

To  address  the  potential  effects  of  climate  change  on  water  yield  in  the  Lower  Athabasca  Regional  Plan  Area 
(LARP),  which  includes  the  downstream  portion  of  the  Athabasca  River  Basin  and  the  Beaver  River  Basin,  the 
HSPF  model,  calibrated  and  validated  for  these  two  basins,  has  been  used  to  simulate  the  hydrologic  effects  of 
forecasted  future  climate  scenarios.  The  1961  to  1990  period  was  used  as  the  climatological  baseline  period  for 
the  effects  assessments. 

As  a  first  step  in  addressing  the  potential  effects  of  climate  change  on  key  ecosystem  variables  in  Alberta, 
Alberta  Environment  has  developed  a  comprehensive  database  (Alberta  Climate  Model)  of  thirteen  climate 
variables  from  the  available  Environment  Canada’s  climate  records  from  1961  to  1990  (Alberta  Environment 
2005).  AENV  provided  Golder  with  the  Alberta  Climate  Model  data  and  average  monthly  changes  in 
temperature  and  precipitation  predicted  by  five  GCMs  for  two  future  periods:  2010  to  2039  (referred  to  as  2020s) 
and  2040  to  2069  (referred  to  as  2050s).  Average  (1961-1990)  monthly  temperature  and  precipitation  values 
are  available  for  grids  of  size  of  approximately  1  km  by  1  km  (0.0083333°)  covering  the  entire  province.  The  five 
climate  scenarios  (2020s  and  2050s)  provided  by  AENV  were: 

h  CCSRNIES_A1  FI  (warmer  and  drier  than  median  conditions); 

■  CGCM2_B23  (cooler  and  drier  than  median  conditions); 

■  HADCM3_A2A  (warmer  and  wetter  than  median  conditions); 

■  HADCM3_B2B  (median  conditions);  and, 

■  NCARPPCM_A1  B  (cooler  and  wetter  than  median  conditions). 

The  average  (2010-2039  and  2040-2069)  monthly  temperature  and  precipitation  values  are  made  available  for 
the  same  grids  as  the  baseline  climate  for  the  entire  province.  The  baseline  (1961-1990),  2020s  and  2050s 
average  monthly  temperature  and  precipitation  values  for  each  of  the  sub-basins  included  in  the  HSPF  model 
were  estimated  as  the  average  of  the  grid  cell  values  within  the  sub-basins.  The  differences  between  the 
baseline  averages  and  the  2020s  averages,  and  the  differences  between  the  baseline  averages  and  the  2050s 
averages  were  then  estimated.  The  average  change  in  mean  annual  precipitation  varies  from  -1%  to  +7%  for  the 
2020s  scenarios  and  from  +2%  to  +13%  for  the  2050s  scenarios.  The  range  of  the  change  in  precipitation  is 
much  wider  on  a  monthly  basis:  -10%  to  +23%  for  the  2020s  scenarios  and  -15%  to  +46%  for  the  2050s 
scenarios.  The  average  increase  in  mean  annual  temperature  varies  from  0.63°C  to  1.52°C  for  the  2020s 
scenarios  and  from  1 ,77°C  to  4.35°C  for  the  2050s  scenarios. 

The  Alberta  Climate  Model  provides  only  the  average  (1961-1990)  monthly  temperature  and  precipitation  values. 
The  HSPF  model  is  a  continuous  simulation  model  that  requires  temperature  and  precipitation  inputs  as  daily 
values.  Using  the  same  data  from  which  the  Alberta  Climate  Model  was  derived,  daily  series  of  temperature 
and  precipitation  were  compiled  from  data  at  seven  index  climate  stations  in  the  LARP  area.  The  observed  data 
were  used  with  the  estimated  differences  for  the  2020s  and  2050s  to  create  daily  series  of  future  climate 
scenarios,  which  were  then  used  with  HSPF  to  develop  forecasted  flow  series  for  the  2020s  and  2050s. 

The  HSPF  model  calibrated  for  the  Athabasca  River  Basin  and  Beaver  River  Basin  was  run  with  the  baseline 
climate  data  and  the  adjusted  future  climate  data.  The  effects  on  flows  were  assessed  using  the  flow  statistics  at 
selected  locations  on  the  main  stem  of  the  Athabasca  River  and  on  the  Beaver  River  at  the  Cold  Lake  Reserve 
for  the  2020s  and  2050s  climate  scenarios. 
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General  conclusions  are  as  follows: 

■  All  the  five  climate  scenarios  for  the  2020s  tend  to  result  in  lower  mean  annual  flows,  with  the  decrease  in 
mean  annual  flow  becoming  more  severe  downstream  along  the  Athabasca  River. 

■  Changes  in  mean  monthly  flows  for  the  2020s  scenarios  tend  to  be  positive  (increase  relative  to  baseline 
values)  during  the  spring  months  and  significantly  negative  (decrease)  during  the  summer  months.  The 
seasonal  differences  in  flows  are  consistent  with  the  relative  effects  of  increased  precipitation  (more  spring 
runoff)  and  increased  temperatures  during  the  summer  months  (higher  summer  evapotranspiration). 

■  The  2050s  climate  scenarios  result  in  changes  in  mean  annual  flow  that  are  generally  more  severe 
compared  to  the  2020s  scenarios,  which  is  consistent  with  the  forecasted  changes  in  temperature  and 
precipitation. 

■  The  range  of  effects  of  future  climate  regimes  on  flows  in  the  Lower  Athabasca  Regional  Plan  Area  (LARP) 
is  summarized  below.  Three  locations  are  used  for  the  summary:  (1)  Athabasca  River  at  Athabasca 
(07BE001 )  to  represent  inflow  to  the  LARP  area,  (2)  Athabasca  River  below  McMurray  (07DA001)  to 
represent  flows  in  the  oil  sands  region  where  the  amount  of  water  withdrawal  is  a  key  consideration  in  the 
planning  process,  and  (3)  Beaver  River  at  Cold  Lake  Reserve  sub-basin  (06AD006)  to  represent  the 
southern  portion  of  the  LARP  area. 

The  summary  results  indicate  that  water  yield  from  the  two  basins  will  generally  decrease,  with  the  effects 
(in  percentage  terms)  more  significant  during  the  month  of  August.  The  effects  on  low  February  flows, 
although  reduced,  tend  to  be  less  (in  percentage  terms)  compared  to  changes  in  August  flows.  Flood  flows 
tend  to  be  higher  under  the  future  climate  scenarios. 

Range  of  Potential  Future  Climate  Effects  on  Flows  in  the  Lower  Athabasca  Region 


Flow  Statistic 

Athabasca  River  at  Athabasca 

Athabasca  River  at  McMurray 

Beaver  River  at  Cold  Lake 

Range  of  Change  in  Flow  (%) 

Range  of  Change  in  Flow  (%) 

Range  of  Change  in  Flow  (%) 

2020s 

Low 

High 

Median 

Low 

High 

Median 

Low 

High 

Median 

Mean  Annual  Flow 

-16% 

-5% 

-9% 

-15% 

-2% 

-6% 

-22% 

0% 

-7% 

Mean  August  Flow 

-19% 

-11% 

-16% 

-21% 

-6% 

-14% 

-30% 

-2% 

-22% 

Mean  February  Flow 

-15% 

-4% 

-7% 

-12% 

-1% 

-4% 

-26% 

-3% 

-10% 

Mean  June  or  April  Flow 

-17% 

7% 

-4% 

-17% 

7% 

-5% 

-4% 

11% 

3% 

10-year  Flood  Flow 

-17% 

2% 

-7% 

-18% 

4% 

-3% 

-22% 

17% 

10% 

2050s 

Low 

High 

Median 

Low 

High 

Median 

Low 

High 

Median 

Mean  Annual  Flow 

-21% 

-4% 

-12% 

-21% 

-4% 

-8% 

-19% 

-2% 

-10% 

Mean  August  Flow 

-34% 

-10% 

-24% 

-31% 

-12% 

-21% 

-27% 

-15% 

-24% 

Mean  February  Flow 

-17% 

0% 

-8% 

-13% 

1% 

-8% 

-19% 

-7% 

-12% 

Mean  June  or  April  Flow 

-12% 

2% 

-5% 

-15% 

4% 

-3% 

-5% 

18% 

5% 

10-year  Flood  Flow 

-11% 

-3% 

-7% 

-15% 

-2% 

-4% 

-11% 

24% 

-3% 

Labels  for  Climate  Scenarios 

CCSRNIES-A1F1  Cl 

CGCM-B23  |C2 

HADCM3-A2A  HI 

HADCM3-B2B  H2 

NCARPCM-A1B  N 
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The  results  of  the  assessment  of  the  effects  of  climate  change  on  flows  in  the  Athabasca  River  and  Beaver  River 
basins  tend  to  be  in  general  agreement  with  recent  studies  carried  out  by  independent  researchers.  Studies  on 
the  effects  of  potential  climate  change  on  flows  in  watersheds  in  Alberta  include  those  by  Kerhoven  and  Gan 
(2005)  on  the  Athabasca  River  Basin  (ARB)  and  Pietroniro  et  al.  (2006)  on  the  South  Saskatchewan  River  Basin 
(SSRB).  Kerkhoven  and  Gan  (2006)  predicted  that  mean  annual  flows  in  the  ARB  would  decrease  by  almost 
25%  by  the  last  30  years  of  the  century.  Pietroniro  and  Toth  (2006)  predict  general  reduction  in  flows  in  the 
SSRB  ranging  from  -13%  to  -4%. 

The  range  of  effects  predicted  by  the  HSPF  model  for  the  LARP  area  is  generally  consistent  with  the  results  of 
similar  studies  in  Alberta.  Depending  on  the  flow  statistic,  the  range  of  effects  can  span  more  than  10%,  which 
reflects  the  assumptions  inherent  in  the  formulation  of  the  climate  scenarios  as  well  as  uncertainties  in  the 
predictions.  The  range  of  hydrologic  effects  predicted  by  the  HSPF  model  for  the  LARP  area  reflects  the  range 
of  conditions  that  basin  planners  should  consider.  Hence,  the  predicted  effects  using  HSPF  and  future  climate 
scenarios  provide  a  rational  basis  for  water  management  and  planning  purposes  in  the  LARP  area. 

The  scope  of  work  for  this  study  included  an  assessment  of  flow  variability  relative  to  effects  due  to  forecasted 
climate  scenarios.  Flow  series  can  be  dominated  by  departures  from  mean  conditions  at  seasonal  to  multi- 
decadal  time  scales.  The  variability  has  been  linked  to  sea  surface  temperature  anomalies  in  the  tropical  and 
extra-tropical  regions  of  the  Pacific  Ocean,  ENSO  and  PDO,  respectively.  Wet  and  dry  cycles  have  significant 
implications  for  prairie  communities  and  economies  and  the  management  of  land  and  water  resources.  The 
calibrated  HSPF  model  calibrated  for  the  Athabasca  River  Basin  was  used  on  the  Firebag  River  sub-basin  to 
generate  a  baseline  (1961-1990)  as  well  as  a  54-year  (1954-2007)  series  of  flows  using  available  data  at  the 
Fort  McMurray  Airport  climate  station.  It  was  assumed  that  the  simulated  55-year  flow  data  series  may  capture 
the  variability  due  to  past  ENSO  events  and  possibly  some  past  PDO  events  to  a  greater  extent  than  the 
standard  30-year  (1961-1990)  series.  The  effect  of  future  climate  scenarios  on  flow  variability  was  assessed  by 
comparing  the  standard  deviation  and  coefficient  of  variation  for  the  two  (1961-1990  and  1954-2007)  annual  flow 
series  under  the  baseline,  2020s  and  2050s  climate  scenarios.  The  calculated  coefficients  of  variation  suggest 
that  there  is  no  difference  (difference  not  statistically  significant)  between  the  baseline  1961-1990  and  1954- 
2007  series,  as  well  between  the  same  two  series  under  the  five  climate  scenarios  for  the  2020s  and  2050s. 
One  reason  for  the  lack  of  difference  between  the  two  series  (30-year  and  54-year)  is  likely  because  the  54-year 
is  still  too  short  and  does  not  include  the  extreme  wet  and  dry  cycles  of  the  past  as  suggested  by  analysis  of 
tree-ring  records. 
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1.0  INTRODUCTION 
1.1  Background 

The  Alberta  Water  for  Life  Strategy  was  developed  to  achieve  safe,  secure  drinking  water,  healthy  aquatic 
ecosystems,  and  reliable  water  supplies  of  suitable  quality  for  a  sustainable  economy.  The  potential  gap 
between  water  demand  and  water  supply  in  the  future,  when  water  demands  are  expected  to  increase  and  water 
supplies  may  decrease,  is  a  key  consideration  in  the  Alberta  Water  for  Life  Strategy.  One  of  the  factors  affecting 
future  water  supply  is  the  nature  of  the  future  climate  regime.  The  effects  of  climate  change  on  water  yield  from 
the  basins  of  Alberta  will  affect  water  uses  and  water  management  in  those  basins.  There  is  a  need  to  assess 
the  potential  effects  so  that  watershed  planners  can  adapt  their  plans  to  take  advantage  of  positive  effects  and 
implement  mitigation  measures  to  minimize  the  negative  effects. 

Alberta  is  embarking  on  a  regional  planning  process  and  the  effect  of  climate  change  on  watershed  responses  is 
a  key  consideration  in  that  process.  One  consideration  of  this  regional  process  is  that  there  is  consistency  in 
how  the  effects  of  climate  change  from  the  various  basins  are  estimated  and  assessed.  Depending  on  the 
hydrologic  response  to  climate  change,  the  effects  on  water  yield  from  a  basin  can  increase  or  decrease  or  there 
may  be  shifts  in  seasonal  patterns.  Hydrologic  models  have  commonly  been  used  to  assess  the  effects  of 
climate  change  on  watershed  hydrology.  One  approach  is  to  use  simulated  outputs  from  GCMs  for  the  next  few 
decades  as  inputs  into  hydrologic  models  to  assess  the  hydrological  responses  of  the  basins  under  future 
modelled  climate  regimes.  As  part  of  a  regionally  specific  initiative  for  the  Lower  Athabasca  Regional  Plan, 
Alberta  Environment  (AENV)  has  identified  the  need  for  a  hydrologic  model  to  assist  it  in  assessing  the  effects  of 
climate  change  on  water  yield  in  the  Lower  Athabasca  Regional  Plan  area  (LARP).  The  intent  is  that  the 
outcome  of  this  regional  initiative  can  be  used  as  a  prototype  to  support  similar  projects  in  other  regions. 

Further  to  a  Request  for  Proposals  and  the  submission  of  proposal  by  Golder  Associates  Ltd.  (Golder),  Alberta 
Environment  (AENV)  contracted  Golder  to  develop  criteria  for  the  selection  of  an  appropriate  hydrologic  model, 
to  develop  an  approach  to  implement  the  selected  model,  to  apply  the  selected  model  to  the  Athabasca  River 
and  Beaver  River  Basins  using  five  climate  change  scenarios,  and  to  assess  the  results  of  the  simulations. 
Phase  I  of  the  project  consisted  of  two  tasks.  The  first  task  was  the  selection  of  an  appropriate  model  based  on 
criteria  that  are  flexible  enough  to  be  transferable  to  the  various  Planning  Regions  of  the  province.  The  follow-up 
task  was  to  develop  a  strategy  to  implement  the  selected  model  to  the  LARP  area,  and  the  Athabasca  River 
Basin  and  the  Beaver  River  Basin,  in  particular.  Phase  II  of  the  project  consisted  of  the  following  tasks: 
calibration  and  validation  of  the  selected  model  on  the  Athabasca  and  Beaver  River  basins;  implementation  of 
the  calibrated  model  using  baseline  (1961-1990)  climate  data  and  GCM-forecasted  climate  data  for  five  climate 
scenarios  for  the  2020s  and  2050s;  and  assessment  of  effects  of  climate  change  on  flow  variability. 

This  report  presents  the  model  selection  criteria,  the  recommended  model  for  implementation,  the  results  of  the 
calibration  and  validation  of  the  selected  model  on  the  Athabasca  River  Basin  and  the  Beaver  River  Basin,  and 
the  results  of  the  assessment  of  future  climate  scenarios  on  flows  in  the  two  basins. 
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2.0  MODEL  SELECTION 
2.1  Model  Selection  Criteria 

The  model  selection  criteria  can  be  broadly  categorized  in  terms  of  (1)  range  and  mathematical  representation  of 
hydrologic  processes  in  the  model  and  the  temporal  scale  of  the  simulations  (hours  to  months),  (2)  applicability 
of  calibrated  model  parameters  or  current  simulation  modules  even  under  modified  climate  regimes,  (3)  flexibility 
of  model  to  simulate  hydrologic  processes  over  a  wide  range  of  modeling  scale,  i.e,  from  large  basins  to  small 
tributary  basins;  (4)  extent  of  climate,  watershed,  vegetation,  and  soil  data  requirements  for  running  the  model, 
including  the  ability  to  use  remote  sensing  data,  (5)  ability  of  input  data  structure  of  model  to  support  the  format 
of  downscaled  climate  change  scenario  data;  (6)  outputs  of  model  and  their  relevance  to  support  Regional 
Planning  efforts  under  the  Land  Use  Framework  and  generally  to  water  and  watershed  management  in  Alberta; 
and,  (7)  currency  or  peer-acceptance  of  model.  A  number  of  specific  criteria  and  sub-criteria  were  then 
developed  from  these  broad  considerations. 

Table  A.1  in  Appendix  A  lists  the  models  and  the  model  selection  criteria  that  were  used  to  select  a  hydrologic 
model  for  hydro-climate  modeling  in  Alberta.  The  criteria  belong  to  the  following  general  categories: 

1.  Criteria  relevant  to  objectives  of  study  -  hydrologic  effects  of  climate  change 

These  include: 

a.  Ability  to  incorporate  key  watershed  features  (physiographic,  hydrographic,  geologic,  etc.) 

b.  Ability  to  simulate  hydrologic  processes  dominant  in  region  or  watershed  under  consideration 

c.  Ability  to  simulate  hydrologic  processes  over  multi-year  periods 

d.  Ability  to  simulate  hydraulic  features  (channel,  floodplain  storage,  lake  storage,  etc.) 

e.  Ability  to  produce  outputs  at  spatial  scale  of  interest  (grid  to  watershed  scale) 

f.  Ability  to  produce  outputs  at  time  scale  of  interest  (hourly  to  annual),  including  floods,  droughts, 
monthly  yield  and  annual  yield 

g.  Ability  of  model  to  use  format  of  downscaled  climate  data  (grid,  basin,  etc.) 

These  criteria  define  the  high-level  functions  of  the  model  that  are  key  to  addressing  the  objective  of  a  study, 
namely,  the  assessment  of  the  effects  of  climate  change  on  water  yield.  Each  criterion  is  further  defined  in  terms 
of  specific  sub-criteria  that  may  be  region,  climate  or  specific  output  dependent.  For  example,  while  hydraulic 
routing  is  a  key  criteria,  the  ability  to  simulate  ice  jams  or  break-up  is  likely  not  significant  for  water  yield  studies. 
On  the  other  hand,  it  could  be  a  significant  criterion  for  assessing  the  effects  of  climate  change  on  flood 
magnitude  and  timing  where  floods  are  governed  by  such  events. 

A  numerical  ranking  from  5  (Critical)  to  0  (No  value)  determines  the  importance  of  the  sub-criteria  to  the  specific 
study  objective  being  considered.  A  numerical  ranking  from  5  (Very  Well)  to  0  (Not  Simulated)  describes  how 
well  the  model  meets  the  sub-criteria.  The  product  of  the  two  scores  provides  an  indication  of  the  ability  of  the 
model  to  address  the  key  study  objectives. 

2.  Criteria  Defining  Data  Requirements  of  Model 
These  include: 

h.  Climate 

i.  Watershed 

j.  Stream 

k.  Geology 

l.  Vegetation 

These  criteria  determine  whether  the  spatial  and  temporal  resolution  of  the  data  required  by  the  model  matches 
that  of  the  available  data.  A  numerical  ranking  from  5  (Available)  to  0  (Not  Available)  specifies  the  resolution  of 
the  data  available,  while  the  actual  data  resolution  required  by  the  model  is  specified  in  terms  of  1  (Yes)  and 
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0  (No).  The  product  of  the  two  scores  provides  an  indication  of  whether  the  data  requirements  of  the  model 
match  what  are  available. 

3.  Criteria  Determining  the  Degree  of  Sophistication  of  the  Processes  Simulated  in  the  Model 

These  include: 

m.  Rainfall-Runoff 

n.  Snowmelt-Runoff 

o.  Evaporation/Evapotranspiration 

p.  Routing  in  Streams 

q.  Infiltration 

r.  Glacier  melt 

While  criteria  1  .a  to  1  .g  define  the  process  that  are  key  the  study  objectives,  criteria  3.m  to  3.r  define  the  method 
used  to  simulate  the  processes.  The  method  may  range  from  the  simple  to  the  very  complex,  and  one  model 
may  have  options  for  the  user  to  select  between  several  of  the  methods.  A  numerical  ranking  from  5  (Critical)  to 
0  (No  value)  determines  the  importance  of  the  various  approaches  to  the  specific  study  objective  being 
considered.  The  degree  of  sophistication  in  the  approach  that  is  important  to  the  study  objective  may  be 
dependent  on  the  region  where  the  model  is  being  applied.  While  sophisticated  snowmelt  routines  would  be 
required  in  areas  where  snowmelt  is  the  major  source  of  water  supply,  the  same  degree  of  sophistication  may 
not  be  necessary  in  areas  where  snowmelt  is  a  minor  contributor  to  runoff.  A  numerical  ranking  from  5  (Very 
Well)  to  0  (Not  Simulated)  describes  how  well  the  model  meets  the  sub-criteria.  The  product  of  the  two  scores 
provides  an  indication  of  the  ability  of  the  model  to  address  the  key  study  objectives  using  approaches  or 
methods  that  are  deemed  to  be  necessary  for  useful  model  outputs. 

4.  Criteria  Defining  Ease  of  Use  by  Modeler 
These  include: 

s.  Model  availability 

t.  Training/Support 

u.  Operating  System 

v.  Model  Description 

w.  Modification  of  Model  by  User 

These  criteria  determine  whether  there  is  adequate  documentation  on  the  model  and  whether  the  model  is 
useable  in  the  modeller’s  work  environment.  A  numerical  ranking  from  5  (Critical)  to  0  (No  value)  determines  the 
importance  of  the  various  requirements  to  the  user.  The  ability  of  the  model  to  meet  these  requirements  is 
specified  in  terms  of  1  (yes)  and  0  (no).  The  product  of  the  two  scores  provides  an  indication  of  whether  the 
user  documentation  requirements  match  the  model  specifications. 

The  next  three  set  of  criteria  define  the  ability  of  the  model  to  address  key  user  requirements  during  the  model 
implementation  stage,  such  as  model  calibration,  model  data  management  and  assessment  of  uncertainty  in 
model  outputs.  These  criteria  are: 

5.  Criteria  for  Model  Calibration  and  Validation 

These  criteria  determine  whether  the  model  has  in-built  modules  or  can  output  simulation  results  that  facilitate 
the  assessment  of  calibration  or  validation  performance.  Sub-criteria  include  graphics  modules  to  display 
observed  and  simulated  stream  flow  series,  generate  statistics  of  observed  and  simulated  series,  and  for 
automatic  optimization  of  calibration  parameters  based  on  the  model  performance  statistics. 

6.  Criteria  for  Data  Management 

These  criteria  determine  whether  the  model  has  user-friendly  input  and  output  data  management  systems.  Sub¬ 
criteria  include  ability  of  model  to  import  input  data  from  standard  spreadsheets  such  as  Excel  and/or  from  GIS 
files  or  remotely  sensed  data  files.  The  ability  to  the  model  to  format  the  simulation  results  for  output  to  similar 
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data  formats  is  also  a  sub-criterion.  The  internal  data  management  capabilities  of  the  model  can  also  be  a  key 
user  requirement. 

7.  Criteria  for  Uncertainty  Analysis 

These  criteria  determine  whether  the  model  has  the  capability  for  uncertainty  analysis  on  simulation  results, 
such  as,  generation  of  random  values  of  calibrated  parameters  from  user-specified  probability  distributions  of 
expected  deviations  of  calibration  parameters  from  their  expected  calibrated  values  or  ability  to  accept  multiple 
synthetically-generated  input  data  series. 

A  numerical  ranking  from  5  (Critical)  to  0  (No  value)  determines  the  importance  of  the  various  requirements  to 
the  user.  The  ability  of  the  model  to  meet  these  requirements  is  specified  in  terms  of  1  (yes)  and  0  (no).  The 
product  of  the  two  scores  provides  an  indication  of  whether  the  model  meets  these  requirements. 

The  sum  of  all  the  scores  for  each  model  is  used  as  an  indication  of  the  appropriateness  of  the  various 
candidate  models  to  meet  the  objectives  of  the  study,  to  have  the  necessary  data  to  run,  as  well  as  to  meet  user 
requirements. 

2.2  Models  under  Consideration 

Tables  A.2  to  A.1 1  in  Appendix  A  describe  the  assessments  of  the  ten  hydrologic  models  that  were  evaluated  for 
the  purposes  of  this  study:  WATFLOOD,  MISBA,  MIKE-SHE,  HSPF,  HEC-HMS,  HYDROTEL,  SSARR,  DHM- 
RS,  VIC  and  CHRM. 

2.2.1  Conceptualizing  Hydrologic  Processes  in  a  Hydrologic  Model 

A  hydrologic  model  transforms  input,  l(t),  to  output,  Q(t),  via  a  set  of  equations  or  transfer  function,  /,  such  that 

Q(t)=/[i(t)]  (i) 

The  solution  of  Equation  (1)  uses  the  basic  principles  of  conservation  of  mass,  momentum,  and  energy,  but  does 
not  generally  account  for  the  precise  flow  pattern  inside  a  control  volume.  Given  that  the  Rainfall-Runoff  (r-r) 
relationship  is  highly  nonlinear,  Equation  (1)  would  be  too  complex  to  characterize  all  the  components  of  the 
rainfall-runoff  transformation  process.  Therefore,  in  modeling  r-r,  modeling  concepts  (model  parameters)  are 
introduced  to  disaggregate  and  simplify  the  complex  process  to  concepts  that  are  manageable,  for  example: 

(1)  Threshold  Concept  of  Infiltration 

Surface  Runoff  r=(i-f)  if  (i-f)  >0  where  /  =  rainfall  intensity  and  f  =  rate  of  infiltration  which  is  function  of  soil 
moisture. 
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(2)  Variable  Source/Contributing  area  Concept 


Dunne  and  Leopold  (1978,  Fig.  9-1 1) 


The  %  of  basin  area  contributing  to  runoff  during  the  course  of  a  storm  varies  with  time.  As  the  water  table 
approaches  the  surface,  the  contributing  area  increases  and  vice  versa,  e.g.,  coupled  effects  of  surface  cover 
and  soil  moisture  (see  diagram  of  Dunne  and  Leopold,  1978). 

2.2.2  Modeling  Scales  and  Approaches 

Mathematical  models  are  constructed  at  some  preferred  scales  for  conceptual  analysis,  for  verification  by 
observations,  and  for  predictive  purposes.  Outputs  from  the  same  models,  but  executed  at  different  spatial 
scales,  can  be  somewhat  different  because  assumptions  made  at  one  scale  may  contradict  the  assumptions 
made  at  another  scale.  Spatial  scales  appropriate  for  hydrologic  models  are  generally  of  the  order  of  10  - 
100  km,  with  time  scales  generally  one  day.  In  contrast,  for  atmospheric  models  like  the  General  Circulation 
Models  (GCMs)  used  for  forecasting  future  climate  regimes,  spatial  scales  are  typically  100’s  km,  and  time 
scales  are  minutes  to  hours.  One  common  approach  to  reconcile  the  different  scales  of  most  hydrologic  models 
and  the  GCMs  is  to  parameterize  those  processes  at  the  micro-scale  (sub-basin)  that  affect  water  behaviour  at 
the  meso-scale  (large  watersheds  or  regions)  as  hydrologic  concepts  at  micro-scale  are  well  established,  but  not 
so  at  mesoscale.  Another  approach  is  to  disaggregate  the  macro-scale  concept  to  allow  for  processes  and  for 
parameter  variations  at  the  meso-scale  of  interest,  but  there  are  few  such  attempts,  e.g.,  Schaake’s 
disaggregation  model  from,  annual  to  monthly  time  step. 

2.2.3  Brief  Descriptions  of  Models 

Some  models  are  more  based  on  parameterization,  others  on  physically-based  processes,  depending  on  the 
types  and  resolution  of  data  available  to  drive  the  model.  The  former  models  are  less  data  intensive  than  the 
latter  models,  but  they  also  tend  to  depend  more  on  parameter  calibrations.  The  models  (Tables  A. 2  to  A.11  in 
Appendix  A),  described  briefly  below,  represent  quite  a  wide  range  of  hydrologic  models  developed  in  the  last 
several  decades. 

WATFLOOD 

WATFLOOD  was  developed  by  Kouwen  (1998)  of  the  University  of  Waterloo.  As  a  fully  distributed  model, 
WATFLOOD  models  a  watershed  in  terms  of  grouped  response  units  (GRUs),  and  is  also  designed  to  use 
gridded  data  sources  such  as  land  cover,  DEM’s,  numerical  weather  prediction  model  output,  and  radar  data. 
Since  it  uses  GRUs  and  land  cover  data,  it  is  possible  to  use  “universal  parameters”  with  some  fine  tuning  or 
calibration.  The  model,  written  in  separate  programming  units,  runs  at  hourly  time  steps  and  can  simulate  up  to 
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100  storm  events.  It  divides  a  watershed  into  a  surface  upper  zone  (saturated,  varying  depth),  an  unsaturated 
zone,  and  a  saturated  lower  zone.  Hydrologic  processes  considered  are  evaporation  loss  using  Hargreaves,  or 
Priestley-Taylor  algorithms,  soil  infiltration  by  Green-Ampt’s  model,  Hortonian  runoff,  Eric  Anderson’s  snow 
model  (temperature  index)  with  cover  based  snow  cover  depletion  curves,  Hartly’s  sediment  model  and  a  crude 
glacier  melt  model.  For  routing,  it  uses  a  center  difference,  kinematic  wave  routing  scheme  with  variable  time 
steps  to  satisfy  the  Courant  criteria.  WATFLOOD  has  been  used  with  grid  sizes  ranging  from  1  to  25  km  and  for 
watershed  areas  from  15  to  1,700,000  km2. 

MISBA 

MISBA  is  a  fully  distributed,  soil-vegetation-atmosphere  transfer  model  (SVAT)  that  considers  interactions 
between  the  soil-biosphere-atmosphere  (ISBA)  of  Me'te'o  France  (Noilhan  and  Planton,  1989),  and  modified 
and  applied  to  the  Athabasca  River  Basin  (ARB)  to  model  its  water  and  energy  fluxes  by  Kerkhoven  and  Gan 
(2006).  It  models  the  energy  and  water  processes  at  the  land  surface  directly  on  the  basis  of  eight  parameters 
at  time  steps  in  seconds:  soil  clay  content,  soil  sand  content,  fraction  of  land  covered  by  surface  water,  depth  of 
soil  column,  heat  capacity  of  vegetation,  and  three  stomatal  resistance  parameters.  Its  input  data  includes  ERA- 
40  historical  re-analysis  dataset  developed  by  the  European  Centre  for  Mid-range  Weather  Forecasts  (ECMWF), 
DEM  data,  landuse  types,  initial  soil  moisture,  initial  soil  ice  content,  initial  deep  soil  temperature,  surface  albedo, 
leaf  area  index,  vegetative  cover  fraction,  vegetative  roughness  length,  and  radiative  fluxes.  It  considers  the 
sub-grid  heterogeneity  of  the  moisture  capacity  of  the  soil  to  follow  the  Xinanjiang  distribution. 

MIKE-SHE 

MIKE-SHE,  developed  by  Abbott  et  al.  (1986),  offers  several  different  approaches  for  hydrologic  simulation, 
ranging  from  simple,  lumped  and  conceptual  approaches  to  advanced,  distributed  and  physically-based 
approaches.  It  considers  precipitation  (rain  or  snow),  evapotranspiration,  including  canopy  interception,  overland 
sheet  flow,  channel  flow,  unsaturated  sub-surface  flow,  and  saturated  groundwater  flow.  If  the  unsaturated  zone 
is  considered,  MIKE-SHE  will  calculate  infiltration,  actual  evapotranspiration  and  recharge;  otherwise, 
groundwater  recharge  must  be  specified.  The  unsaturated  zone  is  a  vertical  soil  profile  model  that  interacts  with 
both  the  overland  flow  and  the  groundwater  model.  Infiltration  can  be  computed  using,  either,  a  simple  2-layer 
root-zone  mass  balance  approach,  a  gravity  flow  model,  or  a  full  Richard’s  equation  model,  all  of  which  require 
specification  of  certain  soil-properties.  MIKE-SHE  comes  with  a  database  of  various  soils  and  typical  crops  for 
different  climatic  regions  of  the  world,  either  through  pseudo-transfer  functions  that  link  suction,  water  content 
and  hydraulic  conductivity,  or  soil  properties  specified  directly.  Groundwater  flow  is  computed  using  a  2D  or  3D 
finite-difference  groundwater  model  similar  to  MODFLOW.  The  river  flow  component  is  a  dynamic,  1  -D  modeling 
tool,  ranging  from  a  simple  Muskingum  routing  to  the  higher  order  dynamic  wave  formulation  of  the  Saint-Venant 
equations,  while  the  overland-flow  component  includes  a  2D  finite  difference  diffusive  wave  approach  using  the 
same  2D  mesh  as  the  groundwater  component.  The  overland  flow  module  is  designed  to  interact  with  the  river, 
the  unsaturated  zone,  and  saturated  groundwater  zone. 

HSPF 

HSPF  uses  a  version  of  the  Stanford  Watershed  Model  (Crawford  and  Linsley,  1966),  which  is  a  deterministic, 
lumped,  conceptual  hydrologic  model  and  so  it  has  modest  data  requirements.  It  distributes  the  incoming  rainfall 
into  canopy  interception,  impervious  areas,  upper  zone  which  will  appear  as  surface  runoff  or  interflow,  and 
infiltration  into  the  lower  zone  storage  divided  into  active  and  inactive  groundwater  storages.  The  three 
conceptual  storages  regulate  soil  moisture  and  groundwater  storages,  while  evapotranspiration  can  extract 
moisture  from  the  interception,  upper,  lower  and  groundwater  storages.  Runoff  from  the  channel  inflow  is  routed 
by  a  hydrologic  routing  technique  that  accounts  for  attenuation  by  the  storage  effect  of  the  channel. 

HEC-HMS 

HEC-HMS  is  a  new  version  of  HEC-1  developed  by  the  US  Army  Corps  of  Engineers’  Hydrologic  Engineering 
Center  (HEC,  1998).  Over  the  years,  it  has  undergone  many  versions.  It  is  also  a  soil-moisture  accounting 
model.  The  latest  version  includes  synthetic  unit  hydrograph  by  Snyder’s  method,  instantaneous  unit 
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hydrograph  by  Clark’s  method,  SCS  curve  number,  exponential  decay  function,  Holtan’s  equation  and  TR-20  to 
compute  infiltration  loss,  streamflow  routing,  Green-Ampt  infiltration  equation,  degree-day  snowmelt  modeling, 
Muskingum-Cunge  flood  routing,  subarea  runoff  computation,  hydrograph  combining.  Parameters  calibration 
can  be  done  using  the  univariate  gradient  method,  and  Nelder  and  Mead  Simplex  method. 


HYDROTEL 


HYDROTEL,  developed  by  Fortin  et  at.  (1986),  is  a  spatially  distributed  hydrological  model  based  on  physical 
processes  specifically  developed  to  facilitate  the  use  of  remote  sensing  and  geographical  information  system 
data.  Various  processes  are  written  in  modules  to  allow  users  having  a  choice  of  algorithms  with  respect  to  the 
data  available  to  drive  the  model,  even  allowing  addition  of  new  algorithms  or  modules.  The  model  can  run  at 
various  time  steps  from  one  to  24  hours,  and  it  consists  of  the  following  components:  meteorological  data, 
accumulation  and  melt  of  the  snow  cover,  potential  evapotranspiration,  vertical  water  budget,  surface  and  sub¬ 
surface  runoff,  and  river  routing. 


SSARR 


Developed  in  1966  by  US  Army  Corps  of  Engineers,  the  SSARR  model  has  been  used  in  many  locations  around 
the  world  including  North  America,  Central  and  South  America,  Europe,  Southeast  Asia,  and  the  Indian 
peninsula  under  a  variety  of  climatic  conditions  and  basin  sizes  ranging  from  5  to  over  400,000  km2.  The 
SSARR  model  is  a  computerised  watershed  model  plus  a  river  and  reservoir  system  model.  It  was  developed  to 
provide  hydraulic  simulations  for  system  analyses  for  planning,  design,  and  operation  of  water  control  works  and 
for  operational  river  forecasting.  The  model  combines  basin  runoff  techniques  producing  surface  and  subsurface 
flow  with  river  and  reservoir  routing  techniques.  Special  optional  routines  include  snow  accumulation  and 
snowmelt,  soil  moisture,  recommended  for  calibration.  Daily  data  from  multiple  stations  include  observed 
precipitation  and  streamflow,  plus  forecasted  precipitation,  and  maximum  air  temperature.  Snow  conditions,  soil 
moisture,  solar  radiation  forecast,  etc.  can  either  be  input  or  allowed  to  be  computed  as  functions  of  other 
parameters.  Printed  output  can  include  all  input  data  plus  forecasted  information  including  surface  and 
subsurface  flow,  percentage  runoff,  snow  water  equivalent,  snowline  elevation,  soil  moisture,  river  stages  and 
discharges,  and  reservoir  inflows,  outflows  and  storage  volumes. 


DPHM-RS 


DPHM-RS  was  developed  by  Biftu  and  Gan  (2001 ,  2004)  for  the  semi-arid  climate  of  the  Canadian  Prairies.  The 
model  can  adequately  account  for  a  river  basin’s  terrain  features  by  sub-dividing  it  into  sub-basins  of  uneven 
shapes  and  sizes  (semi-distributed)  based  on  topographic  information  derived  from  the  digital  terrain  elevation 
(DTED)  data.  It  is  computationally  modest,  designed  to  effectively  assimilate  remotely  sensed  (RS)  data,  and 
has  most  of  its  parameters  determined  through  RS  data  and  measurements.  The  hydrologic  processes  are 
estimated  for  each  land  cover  and  then  aggregated  according  to  the  percentage  of  each  land  cover  present 
within  each  sub-basin.  Evapotranspiration  (ET)  from  each  land  cover  is  estimated  at  three  levels  by  the  two- 
source  model  that  separately  considers  evaporation  from  soil  and  plants.  The  soil  moisture  at  the  top  active  and 
the  transmission  zones  are  estimated  by  a  water  budget  approach,  while  the  groundwater  dynamics  by  the 
topographic  soil  index  obtained  from  DTED.  The  surface  runoff  from  each  sub-basin  is  routed  to  the  channel 
network  by  a  kinematic  wave  response  function,  and  then  routed  to  the  basin  outlet  by  the  Muskingum-Cunge 
model. 

vsc 

VIC-3L  is  a  three-layer  Variable  Infiltration  Capacity  (VIC-3L).  It  is  a  hydrologically  based  land  surface  model 
that  is  classified  under  the  surface  vegetation  atmospheric  transfer  scheme  (SVATS).  VIC-3L  considers  both 
energy  and  water  balances,  and  represents  explicitly  the  effects  of  multiple  vegetation  covers  on  water  and 
energy  budgets.  Input  data  are  such  as  vegetation  covers,  soil  classifications,  and  forcing  data  such  as 
precipitation  and  air  temperature  data.  VIC-3L  also  incorporates  the  representation  of  sub-grid  spatial  variability 
of  precipitation  with  the  representation  of  spatial  variability  of  infiltration  to  simulate  energy  and  water  budgets 
(e.g.,  energy  fluxes,  runoff  and  soil  moisture).  It  includes  both  the  saturation  and  infiltration  excess  runoff 
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processes  in  a  model  grid  cell  with  consideration  of  the  subgrid-scale  soil  heterogeneity,  and  the  frozen  soil 
processes  for  cold  climate  conditions.  Liang  and  Xie  (2001)  improved  VIC-3L’s  (Liang  et  al.,  1996)  runoff 
predictions  by  developing  a  runoff  parameterization  scheme  incorporating  both  Dunne  and  Horton  runoff 
mechanisms. 

CRHM 

CHRM  model  is  a  physically-based  model  primarily  developed  by  the  University  of  Saskatchewan  for  cold  region 
applications  and  muskeg  hydrology  in  northern  Saskatchewan.  Using  physically-based  algorithms,  CRHM 
considers  blowing  snow  transport  by  wind,  snow  interception,  sublimation,  snowmelt,  infiltration  into  frozen  soils, 
hillslope  water  movement  over  permafrost,  actual  evaporation,  and  radiation  exchange  to  complex  surfaces, 
overland  flow,  organic  layer  subsurface  flow,  mineral  interflow,  groundwater  flow,  and  streamflow.  CRHM  has  no 
provision  for  calibration;  parameters  and  model  structure  are  selected  based  on  the  understanding  of  the 
hydrological  system.  It  models  a  watershed  using  the  concept  of  Hydrologic  Response  Unit  (HRU).  More 
information  about  CRHM  is  given  in  Pomeroy  et  al.  (2007). 

2.3  Selected  Model 

Table  A.1  in  Appendix  A  shows  each  model’s  total  score.  The  detailed  evaluation  of  each  model  against  the 
selection  criteria  are  provided  in  Tables  A.2  to  A.11  Appendix  A.  The  model  rankings,  based  on  their  total 
scores  and  suitability  for  this  project,  can  be  categorized  as  follows: 

Top  tier:  MIKE-SHE,  HSPF,  WATFLOOD 

Middle  tier:  DPHM-RS,  VIC,  CRHM 

Bottom  tier:  HYDROTEL,  SSARR,  MISBA,  HEC-HMS 

DPHM-RS,  VIC  and  CRHM  are  sophisticated,  currently  research-oriented,  hydrologic  models  with  modules  that 
tend  to  simulate  hydrologic  processes  using  sophisticated  and  accepted  physically-based  formulations.  Ideally, 
such  models  have  the  potential  to  provide  scientifically  very  defensible  results  on  the  effects  of  climate  change 
on  water  yield.  However,  these  models  are  still  under  development  with  few  or  limited  real-world  applications 
and  uncertain  technical  support.  These  models  may  still  have  “bugs”  or  errors  since  they  are  yet  to  be  tested  in 
wide  usage.  These  models  will  require  significant  effort  and  interaction  with  the  model  developers  over  a  long 
period  of  time  to  implement  province-wide  or  for  multiple  large  watersheds.  Hence,  despite  their  inherent 
strengths  as  physically-based  models  (hence,  their  relatively  high  scores),  they  are  not  recommended  for  this 
project’s  particular  application.  Nevertheless,  a  comparison  between  VIC  and  the  recommended  model  (HSPF) 
is  presented  below  to  highlight  each  other’s  strengths  and  weaknesses.  The  information  can  be  useful  for  a 
future  study  where  the  constraints  of  a  research  model  are  of  lesser  importance. 

The  bottom-tier  models  (HYDROTEL,  SSARR,  MISBA,  HEC-HMS)  are  commercial  or  public-domain  models  that 
have  been  applied  in  real-world  applications,  with  some  of  the  applications  reported  in  either  the  grey-  or  peer- 
reviewed  literature,  and  have  reasonable  levels  of  technical  support.  However,  these  models  are  not  deemed  to 
be  appropriate  for  this  project’s  particular  application  for  the  following  reasons.  MISBA,  a  grid-based  model,  is 
best  suited  for  large  watersheds  and,  once  set  up  for  such  watersheds,  can  be  difficult  to  extract  runoff 
characteristics  at  the  sub-basin  level  because  of  the  lack  of  a  channel  routing  component.  There  is  limited  real- 
world  applications  of  HYDROTEL  reported  in  the  literature,  with  virtually  no  technical  support.  In  addition,  the 
model  codes  and  descriptions  of  the  codes  tend  to  be  mostly  in  French,  hence,  limiting  its  present  application 
within  Alberta  Environment  with  few  francophone  modellers.  SSARR  is  a  conceptual  precipitation-runoff  model 
that  is  very  good  at  flood  forecasting  and  reservoir  routing  simulations,  however,  its  hydrologic  formulations  may 
not  be  sophisticated  enough  to  capture  the  effects  of  climate  change  on  water  yield.  The  snow  melt  routine  in 
HEC-HMS  is  based  on  the  relatively  simple  degree-day  approach.  The  soil  water  balance  and  groundwater 
simulation  modules  are  also  based  on  relatively  simple  algorithms  that  may  not  adequately  capture  the  spatial 
variability  of  yield  in  large  watersheds  with  varied  soil  characteristics.  The  simple  routines  may  not  be 
sophisticated  enough  to  capture  the  effects  of  climate  change  on  the  magnitude  and/or  timing  of  snow  melt 
runoff,  which  is  a  key  component  of  the  yield  from  the  Athabasca  River  basin.  Hence,  these  models 
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(HYDROTEL,  SSARR,  MISBA,  HEC-HMS)  are  not  recommended  unless  models  in  the  top-tier  and  middle-tier 
cannot  be  implemented. 

The  three  top-ranked  models  are  WATFLOOD,  MIKE-SHE  and  HSPF.  MIKE-SHE  is  ranked  first  among  the 
three  because  of  its  physically-based  approaches  for  simulating  hydrologic  processes,  sophisticated  hydraulic 
routing  capabilities,  and  excellent  technical  support  and  documentation.  The  scores  for  WATFLOOD  and  HSPF 
are  not  too  different  from  each  other.  WATFLOOD  is  in  the  top-tier  for  similar  reasons  as  MIKE-SHE,  however, 
WATFLOOD  is  not  yet  as  widely  used  as  MIKE-SHE  and  is  still  a  DOS-based  software.  Both  MIKE-SHE  and 
WATFLOOD  are  proprietary  softwares  with  no  option  for  code  modification  by  the  user.  HSPF  is  the  other 
model  in  the  top-tier.  The  representation  of  hydrologic  processes  in  HSPF  is  not  as  sophisticated  as  those  in 
MIKE-SHE,  however,  they  are  still  fairly  rigorous  and  are  widely  acceptable  and  used.  HSPF  is  a  public  domain 
model.  When  only  selection  criteria  1  ( Criteria  relevant  to  objectives  of  study  -  hydrologic  effects  of  climate 
change)  and  3  ( Criteria  Determining  the  Degree  of  Sophistication  of  the  Processes  Simulated  in  the  Model)  are 
considered,  HSPF  and  WATFLOOD  are  scored  very  close  to  each  other,  with  MIKE-SHE  scoring  higher  than 
either.  When  only  selection  criterion  3  is  considered  on  its  own,  the  three  models  are  ranked  as  follows: 
1:  WATFLOOD,  2:  HSPF,  and  3:  MIKE-SHE. 

Based  on  the  model  selection  criteria  and  evaluations  presented  in  Tables  1  and  2,  full  weighting  of 
considerations  for  the  purposes  of  this  particular  project,  discussions  with  AENV  and  on  the  discussion 
presented  below,  it  is  recommended  that  the  HSPF  model,  although  last  in  the  top-tier,  be  implemented  for  the 
Athabasca  River  and  Beaver  River  basins  in  north-east  Alberta  because  of  its  overall  good  compromise  between 
sophistication,  user-friendly  attributes,  and  wide  use.  MIKE-SHE  and  WATFLOOD  are  the  next  most 
appropriate  models  in  order  of  preference.  These  two  models  are  compared  below.  A  comparison  of  VIC  with 
HSPF  is  also  provided  below  to  assist  AENV  if  it  wishes  to  consider  a  research  model  for  a  future  study  or  a 
pilot/demonstration  study. 

The  rationale  for  recommending  HSPF  as  the  model  for  this  project  is  as  follows: 

Rationale  for  HSPF  Recommendation 

Since  its  initial  development  nearly  twenty  years  ago,  the  HSPF  model  has  been  applied  throughout  North 
America  and  numerous  countries  and  in  various  climatic  regimes  around  the  world.  The  modules  for  simulating 
hydrologic  processes  within  the  model  (e.g.,  snowmelt  runoff,  infiltration,  etc.)  are  reasonably  sophisticated  for 
the  purposes  of  this  project.  Based  on  the  numerous  world-wide  practical  applications,  a  database  of  HSPF 
calibration  parameter  values  have  been  compiled  through  the  joint  sponsorship  of  both  the  U.S.  Environmental 
Protection  Agency  and  the  U.S.  Geological  Survey  as  reference  for  practitioners.  The  availability  of  reference 
calibration  parameters  significantly  improves  the  efficiency  (i.e.,  reduces  time  requirement)  of  model  calibration 
and  validation.  In  addition,  reference  to  the  experience  and  results  of  other  applications  increases  the 
defensibility  of  the  outputs  from  HSPF. 

Experience  with  sophisticated  models  indicates  that  much  of  the  effort  associated  with  hydrologic 
implementation  is  spent  on  data  acquisition  and  management.  Hence,  it  is  essential  that  a  successful 
comprehensive  model  include  a  sound  data  management  component.  The  HSPF  model  software  is  based  on  a 
comprehensive  data  management  system  operating  on  direct  access.  The  HSPF  simulation  modules  draw 
required  input  data  from  time  series  stored  in  its  data  management  system  and  are  capable  of  writing  output  to 
it.  Because  these  data  transfers  require  very  few  instructions  from  the  user,  the  effort  required  with  setting  up 
this  hydrologic  model  for  large  and  complex  watersheds  is  minimized.  In  addition,  the  HSPF  model  includes  the 
following  tools  which  are  not  necessarily  available  as  a  package  in  other  hydrologic  model  software: 

^  Source  code,  executable  version,  user's  guide,  and  technical  support; 

m  A  windows-based  and  independent  interactive  interface  that  is  fully  integrated  with  a  software  product, 
GenScn,  that  provides  the  ability  to  change  an  HSPF  input  sequence  interactively,  run  the  model,  and 
analyze  results  graphically;  and, 
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■  WDMUtil,  a  tool  used  to  manage  and  create  the  watershed  data  management  files  (WDMs)  that  contain  the 
meteorological  data  and  other  time  series  data  used  by  HSPF. 

Under  the  sponsorship  of  U.S.  ERA  and  USGS,  HSPF  continues  to  undergo  refinement  and  enhancement  of  its 
component  simulation  capabilities,  along  with  user  support  and  code  maintenance  activities. 

One  drawback  of  the  HSPF  model  for  application  with  gridded  climate  data  is  that  HSPF  is  implemented  at  a 
sub-basin  level,  that  is,  it  is  not  based  on  grid  cells  as  other  comparable  models  are.  This  implies  that  the 
gridded  climate  scenario  data  will  require  aggregation  at  the  sub-basin  level,  which  will  require  some  effort  on 
the  part  of  the  user  and  potentially  loss  of  input  data  details.  In  addition,  aggregation  may  reduce  some  of  the 
spatial  variability  in  the  climate  data.  This  is  important  if  the  spatial  variability  of  climate  parameters  is  expected 
to  change  under  projected  climate  change  scenarios.  However,  this  potential  drawback  can  be  mitigated  to 
some  extent  by  subdividing  the  watershed  into  sub-basins  that  are  reasonably  small  that  spatially  variability  is 
less  of  an  issue.  The  user  should  be  careful  not  to  sub-divide  more  than  necessary  because  routing  of  the  runoff 
from  a  large  number  of  sub-basins  without  detailed  hydraulic  information  can  be  another  source  of  uncertainty. 

MIKE-SHE  compared  to  WATFLOOD 

MIKE-SHE  and  WATFLOOD  operate  at  a  grid  level,  which  enhances  the  ability  to  import  gridded  climate 
scenario  data.  Both  models  have  good  modules  for  simulating  key  hydrologic  processes.  Both  models  can  use 
data  from  remote  sensing  and  GIS  as  inputs.  However,  MIKE-SHE  has  the  advantage  of  including  both  simple 
and  advanced  hydrologic  process  descriptions  that  can  operate  at  either  a  grid-scale  or  a  sub-watershed  scale. 
Hence,  the  user  can  choose  the  hydrologic  module  component  that  more  appropriately  makes  effective  use  of 
the  available  data,  achieves  the  purpose  of  the  hydrologic  modelling  exercise  and  provides  the  required  end 
results  to  maximize  computational  efficiency.  This  choice  is  not  available  in  WATFLOOD  since  it  operates  at  a 
grid  level  only  and  has  limited  choices  for  hydrologic  and  hydraulic  simulations.  In  addition,  MIKE-SHE  can  be 
easily  linked  to  other  regional  and  local  scale  models  such  as  MIKE  1 1  for  flood  routing,  MIKE-WQ  for  simulation 
of  water  quality,  and  FEFLOW  for  groundwater  simulation.  It  includes  a  complete  user  interface,  sophisticated 
input  and  output  tools,  including  animations,  compared  to  WATFLOOD  that  operates  in  a  DOS  environment  and 
provides  only  text  output.  The  drawback  of  MIKE-SHE  is  the  significant  amount  of  effort  required  to  compile  the 
extensive  data  required  to  run  the  model.  MIKE-SHE  is  also  a  very  sophisticated  model  with  a  steep  learning 
curve.  It  requires  considerable  time  to  become  familiar  with  and  to  be  implemented  by  a  novice  user. 

HSPF  compared  to  VIC 

The  following  is  a  summary  comparison  between  VIC  and  HSPF: 

■  Spatial  Scale:  VIC  is  macro  scale  model  that  operates  at  grid  scale  with  typical  scale  of  application  of  1/8  to 
2  degrees  latitude  by  longitude  (minimum  grid  size  ~  10  km).  The  HSPF  model  operates  at  sub-basin 
scale. 

■  Spatial  distribution  of  surficial  soil  type:  The  VIC  model  assumes  one  surficial  soil  type  for  each  grid. 
For  the  HSPF  model,  the  modeller  can  sub-divide  the  sub-basin  to  account  for  various  soil  types. 

■  Spatial  variability  in  soil  moisture  holding  capacity:  The  VIC  model  accounts  for  spatial  variability  of  soil 
moisture  holding  capacity  at  the  grid  level,  while  the  HSPF  model  accounts  for  spatial  variation  of  infiltration 
at  sub-basin  level. 

■  Base  flow:  The  VIC  model  simulates  baseflow  using  a  non-linear  model.  HSPF  simulates  groundwater  flow 
using  a  linear  or  non-linear  model. 

■  Runoff:  The  VIC  model  has  two  runoff  components  (Surface  and  Baseflow)  while  HSPF  model  accounts  for 
three  runoff  components  (Surface,  Interflow  and  Baseflow). 
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m  Snow  transport  and  sublimation:  In  2004  the  snow  transport  and  sublimation  model  component  developed 
John  Pomeroy  and  his  team  was  incorporated  within  the  VIC  model  for  application  on  two  watersheds  in 
Alaska  and  NWT.  The  HSPF  model  does  not  have  a  component  to  simulate  these  phenomena. 

4 

m  Snow  distribution:  VIC  can  take  the  effect  of  snow  distribution  with  elevation  band  in  the  grid  for  snow  melt 
simulation.  HSPF  has  only  one  parameter  for  updating  land  covered  by  snow  at  each  simulation  time  step. 

4 

m  Snow  Melt  Runoff  algorithm:  VIC  uses  a  complete  energy  balance  at  the  canopy  and  ground  levels  to 

simulate  snow  melt  runoff.  HSPF  model  uses  both  the  temperature  index  and  energy  balance  methods. 

However,  the  energy  balance  method  includes  some  empirical  equations  for  simulation  of  sensible,  latent  f 
and  long-wave  radiation.  ^ 

■  Effect  of  frozen  soils:  VIC  includes  a  detailed  method  to  account  for  frozen  soil  conditions  and  for  ^ 
permafrost.  HSPF  uses  a  conceptual  parameter  to  account  for  frozen  conditions  in  winter. 

■  Flow  Routing:  VIC  model  uses  "separate  and  independent"  routing  models  for  routing  flows  from  grid  to  grid  , 

and  through  the  channel  systems.  The  overland  routing  is  based  on  a  "Unit  Hydrograph"  approach  and  this 
needs  to  be  developed  at  the  implementation  stage  and  provided  as  input.  Channel  routing  is  based  on  the  t 

linearized  Saint-Venant  Equation.  The  HSPF  model’s  overland  flow  routing  is  based  on  Manning's  f 

Equation,  while  channel  routing  is  based  on  storage  routing. 

f 

■  Overall  Model  Comparison:  HSPF  has  the  advantage  in  terms  of  practical  application  for  more  than  20 

years,  having  central  data  base  management,  easy  to  set-up  and  has  GUI  for  pre-  and  post¬ 

processing,  Windows-based  application  and  detailed  user  manuals  and  documents.  The  VIC  model  has  t 
several  sub-modules  that  a  user  needs  to  be  familiar  with  and  compile  individually  to  generate  the 
executable  code,  requires  significant  pre-processing  of  input  data  with  no  central  data  management 
systems  and  a  significant  level  of  effort  for  the  two  separate  routing  models.  r 


Summary  of  Model  Selection  Task  ( 

HSPF  is  considered  to  be  a  practical  model  with  good  technical  support  and  sufficient  scientific  rigour  in  its  f 

representation  of  hydrologic  processes  to  address  the  requirements  of  this  project.  It  is  recommended  for  this  - 

project. 

( 

The  MIKE-SHE  and  WATFLOOD  models  are  good  alternative  models,  however,  based  on  the  full  weighing  of 
considerations,  options  and  suitability  for  the  tasks  for  this  particular  project,  HSPF  was  selected  in  the  final 
consultations.  VIC  is  a  sophisticated  research  model  that  AENV  may  consider  for  a  future  study  or  a  pilot  study. 

AENV  gave  Golder  its  approval  to  proceed  with  the  implementation  of  the  HSPF  model  on  the  Athabasca  and 
Beaver  watersheds  for  the  climate  change  modeling.  AENV  noted  that  it  has  some  concerns  with  the  approach 
that  HSPF  uses  to  simulate  the  snowmelt  process.  AENV  requested  that  Golder  provides  it  with  the  details  and 
results  of  the  calibration  and  validation  of  HSPF  as  soon  as  that  is  completed  so  that  treatment  of  snowmelt  can  f 
be  discussed.  During  the  next  phase  of  the  project,  Golder  confirmed  with  AENV  that  the  energy  method  was 
being  used  in  the  implementation  of  the  HSPF  model  to  the  Athabasca  River  Basin  and  the  Beaver  River  Basin. 
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3.0  MODEL  IMPLEMENTATION  STRATEGY 
3.1  Implementation  Approach 

Following  the  model  selection  task,  the  next  step  in  Phase  I  of  the  project  was  the  development  of  a  strategy  for 
implementing  the  selected  model  on  the  Athabasca  and  Beaver  River  Basins  (henceforth  generically  referred  to 
as  “basin”).  The  approach  should  be  flexible  enough  to  accommodate  the  possibility  of  adjustments  that  may  be 
required  for  implementing  a  different  model  selected  for  another  region. 

Essential  elements  of  the  implementation  approach  include: 

■  Compile  available  model  input  data  for  the  basin,  including  data  on  climate,  soil,  vegetation,  etc.,  for 
calibration  and  validation  purposes. 

■  Compile  flow  data  from  gauged  watersheds  for  calibration  and  validation  purposes. 

■  Assess  changes  in  land  use  patterns  in  gauged  watersheds  during  the  period  of  available  climate  and  flow 
data. 

■  Select  one  period  (5  to  10  years)  of  climate  and  flow  data  that  is  relatively  free  of  trends  and  without  land 
use  changes  for  calibration  purposes  and  another  similar  period  (5  to  10  years)  for  validation  purposes. 
The  calibration  and  validation  periods  should  ideally  include  both  hydrologically  wet  and  dry  years. 

■  Discretize  basin  into  sub-basins  or  grids  (depending  on  model  selected)  that  are  physiographically  different, 
have  different  geologic  characteristics,  experience  different  climate  regimes,  and  have  flow  data  at 
convenient  locations  (outlet  of  sub-basin  or  aggregation  of  grids). 

■  Select  climate  stations  whose  data  would  be  representative  of  the  climatic  regime  within  each  gauged 
watershed  selected  for  calibration  and  validation  purposes. 

■  Calibrate  model  for  each  gauged  watershed  within  each  hydrologic  region  encompassed  by  the  basin  (the 
term  hydrologic  region  is  explained  further  below)  using  statistics  such  as  annual  yield,  monthly  yield,  winter 
flows  and  flood  flows.  Compare  observed  and  simulated  hydrographs  visually. 

■  Validate  model  for  each  gauged  watershed  within  each  hydrologic  region  using  input  from  selected  period 
and  comparing  observed  and  simulated  data  for  the  same  period.  Adjust  calibration  parameters  if 
necessary. 

■  Use  model  calibration  parameters  for  each  gauged  watershed  to  represent  model  parameters  for  each  sub¬ 
basin  encompassing  the  gauged  watershed(s)  used  for  calibration  purposes. 

■  Run  the  model  for  the  entire  basin  using  available  climate  data  at  climate  stations  in  the  basin  and  compare 
with  flows  recorded  at  several  hydrometric  stations  on  the  main  stem  of  the  river.  There  are  several  such 
stations  on  the  Athabasca  from  Jasper  to  Fort  McMurray. 

■  If  the  selected  model  uses  sub-basins  for  simulation  and  routing,  the  baseline  (1961-1990)  climate  data 
may  require  aggregation  for  each  sub-basin  represented  in  the  model  schematic.  If  the  selected  model 
works  on  a  grid  basis,  then  the  discretisation  of  each  sub-basin  should  consider  the  grid  size  of  the  climate 
data  coverage. 

■  Run  the  model  for  the  entire  basin  using  the  baseline  (1961-1990)  climate  data  for  the  basin  and  compare 
with  flows  recorded  at  several  hydrometric  stations  on  the  main  stem  of  the  river  over  the  same  time  period 
(i.e. ,  1961-1990).  Summarize  key  hydrologic  variables  (annual  and  monthly  water  yield  particularly)  at  these 
locations  and  compare  with  those  using  the  climate  station  data. 
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■  AENV  will  provide  five  future  climate  scenarios  derived  from  GCMs  using  statistical  downscaling  techniques 
with  the  Alberta  Climate  Model  (which  represents  baseline  conditions)  for  two  future  30-year  time 
snapshots.  Run  the  selected  and  calibrated  hydrologic  model  on  the  Athabasca  River  and  Beaver  River 
Basins  with  the  five  climate  change  scenarios. 

m  Summarize  key  hydrologic  variables  (annual  and  monthly  water  yield  particularly)  at  selected  locations 
within  the  main  stem  of  the  river  in  the  basin.  Compare  with  those  from  using  the  baseline  Alberta  Climate 
Model  data. 

3.2  Detailed  Implementation  Approach 
3.2.1  Compilation  and  Processing  of  Data 

The  selected  model  will  require  baseline  climate,  physiography  and  hydrology  data  for  several  sub-basins  within 
the  basin  for  calibration  and  validation.  The  first  task  in  implementing  the  selected  model  involves  the  collection 
and  analysis  of  the  relevant  data  required  for  the  set-up  of  the  model.  The  information  would  ideally  be  compiled 
within  a  GIS  framework  to  store  the  data  efficiently,  present  the  information  visually,  and  assist  in  interpreting  the 
results  of  the  modeling  analysis. 

Basin  Topographic,  Soil,  Vegetation  Data  - 

■  Create  a  base  GIS  layer  of  the  basin  from  digital  topographic  and  stream  network  maps. 

■  Sub-divide  the  basin  into  a  number  of  sub-basins.  Large  tributary  sub-basins  will  be  further  sub-divided  to 
reflect  the  topographic  and  climatic  variability  within  such  sub-basins.  At  a  minimum,  each  sub-basin 
should  have  at  least  one  gauged  watershed  in  each  hydrologic  region  encompassed  by  the  basin.  For  the 
Athabasca  River  Basin,  the  Clearwater  River  Sub-basin  will  be  represented  as  a  single  sub-basin,  its  size 
notwithstanding,  because  no  hydrologic  region  has  been  defined  for  the  portion  of  the  sub-basin  lying  in 
Saskatchewan. 

Hydrologic  Regions 

As  part  of  the  hydrologic  regionalization  work  completed  for  AENV  by  Golder  in  2006,  the  Province  of 
Alberta  was  classified  into  hydrologic  regions.  The  regionalization  was  based  on  topography,  climate, 
hydrology,  drainage,  geology,  and  soils.  The  delineation  of  the  hydrologic  regions  was  determined  by 
comparing  the  spatial  patterns  in  physiography,  geology  and  climate  with  the  regions  and  on  the  basis  of 
hydrologic  responses.  The  analysis  of  the  hydrologic  response  patterns  and  physiographic-geologic- 
climatic  spatial  patterns  resulted  in  20  hydrologic  regions  for  the  province  that  represent  a  reasonable 
accounting  of  the  various  factors  influencing  the  hydrologic  response  of  a  watershed.  For  example,  there 
are  seven  hydrologic  regions,  from  mountains  to  foothills  to  prairies,  encompassing  the  Athabasca  River 
Basin  up  to  Lake  Athabasca.  See  Figure  1  in  Appendix  9  of  Alberta  Water  Supply  Assessment  Report 
prepared  by  Golder  for  AENV  in  2007.  The  figure  is  reproduced  as  Figure  C.1  in  Appendix  C  of  this  report. 

■  Carry  out  an  initial  visual  analysis  using  GIS  of  the  spatial  variability  of  the  baseline  (1961-1990)  climate 
data. 

h  Overlay  the  baseline  climate  data  on  the  base  basin  GIS  layer  to  identify  regions  where  significant  spatial 
climate  variability  is  evident.  It  is  expected  that  such  spatial  variability  will  already  have  been  captured  by 
using  the  basin’s  hydrologic  regions.  However,  if  further  refinement  is  required,  then  the  sub-basins  will  be 
further  sub-divided  to  reflect  the  spatial  variability  in  climate. 

eb  IT  IS  ASSUMED  THAT  A  SUB-BASIN  LAYER  HAS  BEEN  CREATED  AT  THIS  STAGE. 

■  Use  the  digital  elevation  model  data  to  develop  elevation  bands  or  ground  slope  maps  for  each  sub-basin. 

The  information  is  used  for  delineating  recharge/discharge  areas,  adjusting  climate  data,  routing  modules, 
etc.,  depending  on  the  hydrologic  model  being  implemented. _ _ 
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■  Create  a  GIS  layer  of  the  basin’s  surficial  geology  and  overlay  on  sub-basin  layer  to  obtain  soil 
characteristics  for  each  sub-basin.  The  information  is  used  for  specifying  infiltration  rates,  soil  moisture 
capacity,  etc.,  depending  on  the  hydrologic  model  being  implemented. 

■  Create  a  GIS  layer  of  the  basin’s  vegetation  cover  and  overlay  on  sub-basin  layer  to  obtain  vegetation 
characteristics  for  each  sub-basin.  The  information  is  used  for  specifying  interception  capacity,  estimation 
of  evapotranspiration,  etc.,  depending  on  the  hydrologic  model  being  implemented. 

Climate  Data  -  A  number  of  climate  monitoring  stations  are  in  operation  in  the  basin.  They  include  real-time  and 
daily  recording  stations  operated  by  Environment  Canada,  and  seasonal  stations  operated  by  AENV  to  support 
forest-fire  protection  activities.  A  daily  time  series  of  each  climate  variable  required  by  the  selected  model  should 
be  developed  for  each  station  located  within  or  close  to  gauged  watersheds  that  could  potentially  be  used  in  the 
calibration  or  validation  of  the  selected  model. 

Hydrometric  Data  -  Water  Survey  of  Canada  (WSC)  has  operated  or  continues  to  operate  streamflow 
(hydrometric)  monitoring  stations  on  a  number  of  streams  within  the  basin.  The  gauged  watersheds  cover  a 
range  of  drainage  areas,  elevations,  slopes,  vegetation  characteristics,  etc.  The  data  series  are  not  concurrent, 
but  there  is  good  overlap  of  data  series  from  the  1970s  to  date.  Time  series  of  available  data  on  land  use 
changes  and  other  developments  can  be  compared  to  periods  of  flow  data  to  identify  the  time  periods  that  can 
be  used  for  the  calibration  and  validation  of  the  selected  model. 

Land  Use  Changes  -  Areas  within  the  basin  may  have  over  time  undergone  changes  due  to  natural  resource 
management  activities  (agriculture,  logging,  diversion,  etc.)  and  infrastructure  developments  (roads,  culverts, 
bridges).  For  the  purpose  of  this  study,  an  overview  of  temporal  changes  in  land  use/cover  in  the  sub-basins 
selected  for  calibration  and  validation  would  be  conducted  to  identify  sub-basins  that  are  in  essentially  natural 
conditions. 

Water  Level  in  Lesser  Slave  Lake  -  Lesser  Slave  Lake  is  a  significant  waterbody  in  the  Athabasca  River  basin 
and  Lesser  Slave  River  is  a  major  tributary  to  the  Athabasca  River.  Water  levels  for  Lesser  Slave  Lake  (LSL) 
are  available  from  1979  onwards.  AENV  has  developed  a  “naturalized”  series  of  flow  and  water  level  data,  from 
1916  to  1999,  for  a  location  at  the  outlet  of  LSL.  The  naturalization  was  based  on  the  calibration  and 
implementation  of  a  SSARR  model  for  the  catchment  area  upstream  of  the  lake  outlet.  AENV  also  developed  a 
corresponding  series  of  regulated  flows  and  water  levels  from  1916  to  1999,  based  on  the  known  characteristics 
of  the  weir.  This  is  a  synthetic  series,  assuming  the  presence  of  a  weir  since  1 91 6.  AENV  has  concluded  that  the 
lake  water  level  data  corresponding  to  the  post-weir  recorded  and  synthetic  regulated  series  are  “very  nearly 
coincident”  and  that  the  difference  between  simulated  regulated  levels  and  observed  levels  since  weir 
construction  is  “not  significant”.  The  hydraulic  characteristics  of  the  weir  on  LSL  should  be  incorporated  in  the 
selected  model  to  account  for  any  effects  of  climate  change  on  water  levels  in  the  lake  and  the  consequent 
effects  on  flows  into  the  Athabasca  River. 

River  Cross-Section  Data  -  The  hydraulic  modelling  domain  is  expected  to  include  Lesser  Slave  Lake,  the 
reach  of  the  Lesser  Slave  River  from  the  LSL’s  outlet  to  the  confluence  with  Athabasca  River,  Athabasca  River 
downstream  of  Fort  McMurray  and  up  to,  but  excluding,  Lake  Athabasca,  and  the  reaches  of  other  major 
contributing  tributaries,  contingent  on  the  availability  of  cross-section  data.  Cross-sectional  data  can  be  obtained 
from  the  rating  curves  developed  by  WSC  at  its  gauging  stations. 


3.2.2  Calibration  and  Validation  of  Model 


■  Select  at  least  one  gauged  watershed  with  adequate  flow  records  within  each  hydrologic  region  of  the 
basin.  Select  the  climate  station(s)  to  represent  each  hydrologic  region.  Adjust  temperature  and 
precipitation  data  for  elevation  if  the  climate  station  is  at  an  elevation  significantly  different  from  most  of  the 
sub-basin. 


If  the  model  uses  sub-basins  as  simulation  units,  the  calibration  and  validation  of  the  model  can  proceed  for 
each  selected  gauged  watershed. 
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If  the  model  uses  grids  as  simulation  units,  discretize  each  sub-basin  selected  for  calibration  and  validation 
purposes  into  grid  cells  based  on  homogeneity  in  current  land  use  or  cover;  ground  slope  or  elevation;  soil 
groups;  channel  slope;  etc.  The  digital  elevation,  soil,  land  use,  and  hydrography  data  for  the  basin  will  be 
used  for  the  discretization.  The  extent  of  discretization  or  size  of  the  grid  cells  will  be  a  balance  between  a 
fine  level  of  detail  that  makes  full  use  of  the  physically  based  equations  within  the  selected  model  and  a 
coarser  resolution  that  recognizes  the  limitations  in  data  availability.  The  discretization  will  take  into  account 
land  use  changes  over  time.  Good  representation  of  the  landscape  should  be  achieved  without  excessive 
discretization.  The  compromise  between  discretization  detail  and  data  resolution  implies  that  the  model  will 
require  calibration.  Model  parameters  will  be  adjusted  based  on  model  simulations  and  observed  flow  or 
water  level  data  at  gauged  locations. 

■  Channel  routing  of  runoff  will  be  based  on  cross-section  data,  if  they  are  available,  or  map-derived 
characteristics  (e.g.  slope)  or  area-based  estimates  (i.e.  channel  width  and  depth)  or  generic  values,  if  they 
are  not  available. 

■  Select  two  time  periods  (5  to  10  years  each)  for  which  climate,  hydrology,  and  land  use  data  are  available 
from  gauged  watersheds  under  natural  conditions  (i.e.  no  developments). 

■  Run  the  selected  model  (continuous  simulation)  using  one  of  the  time  periods  to  calibrate  the  model 
(continuous  simulation)  and  the  other  to  validate  the  model.  Adjust  model  parameters  until  calibration  and 
validation  are  deemed  reasonable  based  on  objective  statistical  criteria. 

■  Alternatively,  select  one  gauged  sub-basin  with  a  given  hydrologic  region  for  calibration  and  another  nearby 
gauged  sub-basin  for  validation.  The  climate  data  series  used  for  calibration  and  validation  do  not  need  to 
be  similar  in  this  approach. 

s  The  model  calibration  parameters  for  each  gauged  watershed  can  be  used  for  all  sub-basins  (gauged  and 
ungauged)  in  the  basin,  which  are  similar  in  topographic,  climatic,  soil  and  vegetation  characteristics. 

3.2.3  Running  Climate  Scenarios  with  Calibrated  Model 

The  implementation  steps  following  calibration  and  validation  of  the  model  are  as  listed  above  and  repeated 

hereunder: 

m  Run  the  model  for  the  entire  basin  using  available  climate  data  at  climate  stations  in  the  basin  and  compare 
with  flows  recorded  at  several  hydrometric  stations  on  the  main  stem  of  the  river.  There  are  several  such 
stations  on  the  Athabasca  River  from  Jasper  to  Fort  McMurray. 

h  If  the  selected  model  uses  sub-basins  for  simulation  and  routing,  the  baseline  (1961-1990)  climate  data 
may  require  aggregation  for  each  sub-basin  represented  in  the  model  schematic.  If  the  selected  model 
works  on  a  grid  basis,  then  the  discretization  of  each  sub-basin  should  consider  the  grid  size  of  the  climate 
data  coverage. 

a  Run  the  model  for  the  entire  basin  using  the  baseline  (1961-1990)  climate  data  for  the  basin  and  compare 
with  flows  recorded  at  several  hydrometric  stations  on  the  main  stem  of  the  river  over  the  same  time  period 
(i.e.,  1961-1990).  Summarize  key  hydrologic  variables  (annual  and  monthly  water  yield  particularly)  at  these 
locations  and  compare  with  those  using  the  climate  station  data. 

a  AENV  will  provide  five  future  climate  scenarios  derived  from  GCMs  using  statistical  downscaling  techniques 
with  the  Alberta  Climate  Model  (which  represents  baseline  conditions)  for  two  future  30-year  time 
snapshots.  Run  the  selected  and  calibrated  hydrologic  model  on  the  Athabasca  River  and  Beaver  River 
Basins  with  the  five  climate  change  scenarios. 

a  Summarize  key  hydrologic  variables  (annual  and  monthly  water  yield  particularly)  at  selected  locations 
within  the  main  stem  of  the  river  in  the  basin.  Compare  with  those  from  using  the  baseline  Alberta  Climate 
Model  data. 
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HYDRO-CLIMATE  MODEL  SELECTION  AND  APPLICATION  ON 
THE  ATHABASCA  AND  BEAVER  RIVER  BASINS 


4.0  CALIBRATION  OF  HSPF  MODEL  FOR  THE  LOWER  ATHABASCA 
REGIONAL  PLAN  AREA 

4.1  Introduction 

The  HSPF  model  was  selected  to  assess  the  hydrologic  responses  of  the  Athabasca  River  Basin  and  the 
Beaver  River  Basin  (comprising  the  Lower  Athabasca  Regional  Plan  Area)  under  future  modelled  climate 
regimes  forecasted  by  several  Global  Climate  Models  (GCMs)  for  the  next  few  decades.  This  section  describes 
the  calibration  and  validation  of  the  HSPF  model  for  the  two  basins. 

4.2  Athabasca  River  and  Beaver  River  Basins 

The  Athabasca  River  starts  in  the  Rocky  Mountains  near  Mount  Columbia  (elevation  3,747  metres)  and  flows 
northeast  for  1,400  kilometres  until  it  empties  into  Lake  Athabasca  (elevation  208  metres) 
(http://scienceoutreach.ab.ca/athabasca  reqion.htm.  Accessed  20  June  2009.).  Flows  from  the  basin  eventually 
make  their  way  to  the  Arctic  Ocean.  The  river  drains  an  area  of  approximately  138,000  km2.  The  river  flows 
past  the  urban  centres  of  Jasper,  Hinton,  Whitecourt,  Athabasca  and  Fort  McMurray  prior  to  emptying  into 
Lake  Athabasca.  The  Athabasca  River  Basin  includes  the  McLeod,  Pembina  and  Clearwater  rivers. 

As  a  major  river  system,  the  Athabasca  River  is  influenced  by  a  variety  of  climate,  terrain  and  landscape 
characteristics  found  within  its  basin  (http://www.ramp-alberta.ora/river.aspx.  Accessed  20  June  2009.).  The 
seasonality  of  climatic  conditions  is  a  major  influence  affecting  river  flow  conditions.  The  climate  includes  cold 
winters,  when  most  of  the  seasonal  precipitation  falls  as  snow.  Cold  winters  are  typically  followed  by  warm 
summers,  when  snow  and  glacial  melt  waters  from  the  river’s  headwaters  combine  with  runoff  from  localized 
snowmelt  and  rainfall  events  throughout  the  basin.  As  the  river  flows  toward  Lake  Athabasca,  water  is 
contributed  to  the  river  from  individual  sub-basins. 

The  Athabasca  River  Basin  encompasses  the  following  ecoregions  (natural  regions): 

■  Rocky  Mountain  (Alpine,  Subalpine,  Montain); 

■  Boreal  Foothills; 

■  Boreal  Mixed  Wood  (dry  and  wet); 

■  Boreal  Uplands; 

■  Boreal  Lowlands;  and, 

■  Canadian  Shield  (Athabasca  Plains). 

The  Alberta  portion  of  the  Beaver  River  Basin  is  located  in  the  east-central  part  of  the  province,  near  the  town  of 
Bonnyville.  This  basin  lies  entirely  within  the  Boreal  Mixedwood  Ecoregion  (University  of  Alberta  2005).  The 
Beaver  River  starts  near  the  town  of  Lac  La  Biche  as  the  outflow  from  Beaver  Lake.  The  basin  is  characterized 
by  undulating  to  moderately  rolling  topography,  elevations  ranging  from  about  500  m  to  750  m  above  sea  level 
(AENV  2006)  and  an  abundance  of  lakes.  A  significant  portion  of  the  low  lying  areas  of  the  Beaver  River  Basin 
is  classified  as  non-contributing  because  these  areas  do  not  contribute  stream  flows  during  a  storm  with  a  return 
period  less  than  or  equal  to  10-year.  The  Sand  River  is  the  largest  tributary  to  the  Beaver  River. 

4.2.1  HSPF  Model  Setup  for  the  Athabasca  River  Basin  and  Beaver  River  Basin 

The  continuous  (dynamic)  Hydrologic  Simulation  Program  Fortran  (HSPF)  model  developed  by  the  U.S. 
Environmental  Protection  Agency  (U.S.  EPA)  was  selected  for  the  Athabasca  River  and  Beaver  River  hydro¬ 
climate  modeling  Study.  The  HSPF  model  simulates  stream  flows  as  the  sum  of  three  components:  surface  flow, 
interflow  and  groundwater.  The  relative  magnitude  of  each  component  depends  on  land  use,  soils  and 
vegetation  cover.  The  model  user  can  specify  specific  parameters  for  various  land  use  types  to  represent  the 
physical  processes  in  a  basin. 
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HYDRO-CLIMATE  MODEL  SELECTION  AND  APPLICATION  ON 
THE  ATHABASCA  AND  BEAVER  RIVER  BASINS 


4.2. 1. 1  Characteristics  of  the  Athabasca  River  Basin  and  Beaver  River  Basin 

The  surficial  geology  of  a  basin  is  a  significant  factor  in  the  hydrologic  response  of  the  basin.  HSPF  models  the 
response  of  a  basin  to  precipitation  based  on  the  soil  characteristics  (infiltration  rate,  storage  capacity,  interflow, 
etc.)  of  the  various  land  types  in  the  basin.  In  the  calibration  of  the  HSPF  model  for  the  Athabasca  River  Basin 
and  Beaver  River  Basin,  the  discretization  of  the  sub-basins  into  land  types  was  based  on  surficial  geology.  The 
surficial  geology  characteristics  considered  in  the  set-up  of  the  HSPF  model  for  the  Athabasca  River  Basin  and 
Beaver  River  Basin  are  shown  in  Figure  B.1  in  Appendix  B.  The  land  types  in  the  basin  were  based  on  the 
following  nine  major  surficial  geology  classifications: 

m  Well  Drained  Sand  and  Rapidly  Drained  Sand; 

■  Well  Drained  Till  and  Rapidly  Drained  Till; 

■  Well  Drained  Clay  Loam; 

■  Organic  Soil; 

m  Poorly  Drained  Sand  (Lowland  Glaciolacustrine); 

■  Poorly  Drained  Till  (Lowland  Glaciofluvial); 

■  Poorly  Drained  Clay  Loam  (Lowland  Glacial); 

■  Impervious/Fractured  Rock;  and, 

h  Impervious/Glacier. 

The  headwaters  of  the  Athabasca  River  Basin  are  covered  by  glaciers.  Glaciated  areas  were  modelled  as  an 
impervious  land  type,  but  with  discharges  from  the  glaciated  areas  when  temperatures  rise  above  0°C.  A 
significant  portion  of  the  upper  areas  of  the  Athabasca  River  Basin  is  also  classified  as  impervious.  However, 
most  of  these  areas  are  fractured,  thus  increasing  the  travel  time  of  runoff  compared  to  strictly  impervious  areas. 
Such  areas  were  treated  as  “pervious”  in  the  HSPF  model,  but  with  appropriate  adjustments  of  the  pervious  land 
type  parameters. 

The  characteristics  of  the  vegetation  cover  in  a  basin  also  play  a  significant  role  in  the  basin’s  hydrologic 
response  through  interception  of  precipitation,  evapotranspiration  of  intercepted  precipitation  and  water  stored  in 
the  soil  layer,  shading  of  solar  radiation  (which  affects  snow  melt  rate),  etc.  The  vegetation  characteristics  in  the 
Athabasca  River  Basin  and  Beaver  River  Basin  are  shown  in  Figure  B.2  provided  in  Appendix  B.  The  sub¬ 
basins  in  the  Athabasca  River  Basin  and  Beaver  River  Basin  were  not  explicitly  divided  into  vegetation  types  (as 
was  done  for  land  types  based  on  surficial  geology),  however,  the  characteristics  were  considered  in  the 
specification  of  parameters  controlling  interception,  evapotranspiration  and  shading.  A  sub-division  of  land  types 
based  on  vegetation  can  be  considered  for  future  applications  of  the  model  when  the  effects  of  changes  in 
vegetation  as  a  result  of  climate  change  need  to  be  modelled  can  be  predicted  and  the  hydrologic  effects  need 
to  be  assessed. 

4.2.1. 2  Representation  of  Basins  in  the  HSPF  Model 

The  Athabasca  River  Basin  was  sub-divided  into  75  sub-basins,  based  on  the  drainage  network  and  the 
locations  of  the  WSC  stations  on  gauged  sub-basins  selected  for  calibration  and  validation  of  the  model.  The  75 
sub-basins  were  further  sub-divided  based  on  surficial  geology.  The  locations  of  the  hydrometric  stations  for  the 
selected  gauged  sub-basins  are  shown  in  Figure  C.2  in  Appendix  C.  The  connectivity  between  the  sub-basins 
and  the  characteristics  of  each  sub-basin  were  represented  or  conceptualized  as  shown  in  the  schematic  of  the 
HSPF  model  (Figure  D.1)  in  Appendix  D.  The  schematic  shows  the  “HSPF  model  flow  diagram”,  which  is 
basically  the  sequence  in  which  the  sub-basins  are  simulated  by  HSPF  and  the  connectivity  of  each  sub-basin  to 
downstream  receiving  reach. 

Similarly,  the  Beaver  River  Basin  was  sub-divided  into  three  sub-basins  as  shown  in  Figure  C.2  in  Appendix  C. 
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4.2.1. 3  Model  Data  Input 

Data  used  to  run  the  model  were  as  follows: 

■  Temperature  and  precipitation  data  for  the  calibration  of  the  model  were  from  climate  stations  closest 
(specific  selection  based  on  availability  of  concurrent  climate  and  flow  data)  to  the  gauged  sub-basins.  See 
Figure  C-3  in  Appendix  C.  To  account  for  the  spatial  variability  of  precipitation,  precipitation  data  from 
different  stations  were  assigned  to  the  different  Hydrologic  Regions  encompassed  by  the  Athabasca- 
Beaver  River  Basins.  The  Hydrologic  Regions  of  Alberta  were  developed  by  Golder  (2006)  as  part  of  a 
project  for  AENV.  See  Figure  C.1  in  Appendix  C  for  the  hydrologic  regions  of  the  Lower  Athabasca 
Regional  Plan  Area.  The  Hydrologic  Regions  represent  areas  within  which  the  climate,  geology  and 
hydrologic  responses  are  more  or  less  homogeneous,  but  different  from  the  adjacent  Hydrologic  Region. 
The  Hydrologic  Regions  provide  a  rational  basis  for  assigning  climate  stations  to  the  sub-basins  of  the 
Athabasca  River  Basin  and  Beaver  River  Basin.  Table  4.1  shows  the  climate  stations  selected  for 
precipitation  and  temperature  data  for  the  sub-basins  within  the  Hydrologic  Regions  encompassing  the 
Athabasca  and  Beaver  River  Basins. 


Table  4.1  Climate  Stations  Selected  for  Sub-Basins  within  Hydrologic  Regions  of  Athabasca  and  Beaver 
River  Basins 


Hydrologic  Region  (HR) 

Climate  Station(s)  Selected  for 
Precipitation  and  Temperature  Data  for 
Sub-Basins  within  HR 

Climate  Station(s)  Selected  for  Wind 
Speed  and  Dew  Point  Temperature 
Data  for  Sub-Basins  within  HR 

3 

Jasper  &  Jasper  Warden 

Edson 

4 

Jasper  East  Gate 

Edson 

10.1  (south  of  Athabasca  River) 

Edson 

Edson 

10.2  (north  of  Athabasca  River) 

Slave  Lake 

Edmonton  International  A.  for  Wind  Speed 
Edson  for  Dew  Point  Temperature 

8 

Campsite 

Edmonton  International  A. 

2C 

Lac  La  Biche 

Fort  McMurray 

9A 

Fort  McMurray 

Fort  McMurray 

2D 

Cold  Lake 

Edmonton  International  A. 

■  Wind  speed  and  dew  point  temperature  data  were  from  the  Edmonton  International  Airport  station  or  Edson 
depending  on  the  Hydrologic  Region  as  shown  in  Table  4.1,  while  solar  radiation  data  were  from  the 
Edmonton  Stony  Plain  station.  The  solar  radiation  data  at  the  Edmonton  Stony  Plain  Station  were  assumed 
to  be  representative  of  the  entire  Athabasca-Beaver  River  basins  because  solar  radiation  tends  to  be 
generally  less  spatially  variable  than  most  other  climatic  variables. 

■  Evapotranspiration  and  lake  evaporation  data  were  derived  using  the  Morton  Model,  with  air  temperature, 
dew  point  temperature,  precipitation  and  solar  radiation  used  as  input  data. 

■  Channel  cross-sections  used  in  generating  depth-area-volume-flow  tables  were  estimated  from  data  used 
by  WSC  to  develop  rating  curves  at  the  hydrometric  gauging  stations. 

■  Recorded  stream  flows  at  several  hydrometric  stations  in  the  basin  were  used  for  model  calibration  and 
validation. 
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4.2. 1.4  HSPF  Model  Calibration  and  Validation 

The  approach  to  calibrate  the  HSPF  model  was  modified  slightly  from  that  presented  in  Section  3.2.2.  Instead  of 
calibrating  and  validating  the  model  using  gauged  sub-basins  within  each  Hydrologic  Region,  gauged  sub-basins 
that  are  dominated  by  one  surficial  geology  type  were  selected  and  the  model  parameters  for  this  surficial 
geology  type  “fixed”.  The  process  was  repeated  for  the  other  surficial  geology  types.  The  purpose  of  this 
approach  was  to  establish  “regional  model  parameters”  that  can  be  transferred  to  other  sub-basins  in  the  region 
that  have  similar  surficial  geology.  The  effects  of  land  cover  or  vegetation  types  were  indirectly  considered  in  the 
modelling  through  calibration  of  the  parameters  that  simulate  their  influences.  The  vegetation  types  were 
considered  through  the  following  HSPF  model  parameters:  FOREST  (fraction  of  the  pervious  land  segment  that 
is  covered  by  forest),  CEPSC  (interception  storage  capacity),  and  LZETP  (lower  zone  evapotranspiration 
parameter).  Table  4.2  shows  the  sub-basins  used  for  calibration  and  validation,  as  well  as  location,  drainage 
area,  land  type  and  other  information  on  the  sub-basins. 
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HYDRO-CLIMATE  MODEL  SELECTION  AND  APPLICATION  ON 
THE  ATHABASCA  AND  BEAVER  RIVER  BASINS 


At  least  three  years  of  meteorological  and  hydrologic  data  were  used  to  calibrate  the  HSPF  model.  Calibration 
began  with  an  initial  estimate  of  the  model  parameters  based  the  lower  and  upper  limits  of  each  model 
parameter  as  described  in  the  HSPF  manual.  Then,  the  simulated  monthly  and  annual  runoff  volumes  were 
compared  to  the  observed  volumes  at  hydrologic  gauging  stations.  Appropriate  parameters  were  adjusted  until 
the  simulated  monthly  and  annual  volumes  were  acceptably  close  to  the  observed  values.  The  outcome  of  the 
calibration  and  validation  of  the  HSPF  model  was  deemed  to  be  good  to  reasonable  to  poor  based  on  the 
statistical  criteria  given  in  Section  4.2.1 .5.  Individual  event  volumes  were  then  calibrated. 

Once  streamflow  volumes  were  calibrated,  flow  hydrographs  were  calibrated  using  both  interflow  and  channel 
routing  parameters.  The  shapes  of  event  hydrographs,  and  to  some  extent  the  peak  flows,  were  calibrated  by 
changing  the  interflow  parameters  and  the  appropriate  stage-storage-discharge  relationships.  A  combination  of 
manual  and  automatic  (using  PEST)  calibration  were  used  to  derive  the  model’s  calibration  parameters.  The 
calibration  parameters  for  each  of  the  nine  land  types  are  shown  in  Table  4.3.  Values  in  parentheses  in  Table 
4.3  identify  instances  when  slight  changes  were  made  in  the  parameters  of  similar  land  types  but  located  in 
different  parts  of  the  basin. 

The  specific  model  calibration  approach  for  the  Athabasca-Beaver  River  Basins  is  summarized  as  follows: 

■  Model  parameters  for  the  Well/Rapidly  Drained  Till  land  type  were  calibrated  using  the  recorded  stream 
flows  at  Lovett  River  near  the  Mouth  (Environment  Canada  Hydrometric  Station  07BA003).  The  entire 
Lovett  River  sub-basin  is  covered  with  the  Well/Rapidly  Drained  Till  land  type  (see  Figure  4.1).  The 
precipitation  data  used  were  from  the  Lovett  Lookout  station  and  missing  winter  precipitation  data  were 
filled  using  data  from  the  Edson  climate  station. 

■  Model  parameters  for  the  Well  Drained  Clay  Loam  land  type  were  calibrated  using  the  recorded  stream 
flows  at  Erith  River  below  Hanlan  Creek  (Environment  Canada  Hydrometric  Station  07AF907).  The  surficial 
geology  of  the  Erith  River  sub-basin  is  approximately  60%  Well/Rapidly  Drained  Till  and  40%  Well  Drained 
Clay  Loam  (see  Figure  4.1).  During  the  calibration  process,  the  model  parameters  for  the  Well/Rapidly 
Drained  Till  land  type  as  determined  during  the  calibration  on  the  Lovett  River  sub-basin  were  transferred  to 
the  Erith  River  sub-basin.  The  precipitation  data  used  were  from  the  Lovett  Lookout  station  and  missing 
winter  precipitation  data  were  filled  using  data  from  the  Edson  climate  station. 

■  Model  parameters  for  the  Organic  land  type  were  calibrated  using  the  recorded  stream  flows  at  Mackay 
River  above  Dunkirk  River  (Environment  Canada  Hydrometric  Station  07DB005).  The  entire  Mackay  River 
above  Dunkirk  River  watershed  was  assumed  to  be  covered  by  the  Organic  land  type  (see  Figure  4.2).  The 
precipitation  data  used  were  from  Livock  Lookout  station  and  the  missing  winter  precipitation  data  were 
filled  using  data  from  the  Fort  McMurray  Airport  climate  station. 

■  Model  parameters  for  the  Glacier  land  type  were  calibrated  using  the  recorded  stream  flows  at  Sunwapta 
River  at  Athabasca  Glacier  (Environment  Canada  Hydrometric  Station  07AA007).  The  entire  Sunwapta 
River  sub-basin  was  assumed  to  be  covered  by  glaciers  (see  Figure  4.1).  The  precipitation  data  used  were 
from  Jasper  &  Jasper  Warden  climate  stations. 

■  Model  parameters  for  the  Impervious/Fractured  Rock  land  type  were  calibrated  using  the  recorded  stream 
flows  at  Miette  River  near  Jasper  (Environment  Canada  Hydrometric  Station  07AA001).  The  surficial 
geology  of  the  Miette  River  sub-basin  is  approximately  31%  Well/Rapidly  Drained  Till,  20%  Glacier,  and 
49%  Impervious/Fractured  Rock  (see  Figure  4.1).  During  the  calibration  process,  the  model  parameters  for 
the  Well/Rapidly  Drained  Till  and  Glacier  land  types  were  transferred  from  those  obtained  during  the 
calibration  of  the  Lovett  River  and  Sunwapta  River  sub-basins,  respectively.  The  precipitation  data  used 
were  from  Jasper  &  Jasper  Warden  climate  stations. 
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HYDRO-CLIMATE  MODEL  SELECTION  AND  APPLICATION  ON  THE  ATHABASCA  AND  BEAVER  RIVER  BASINS 


Table  4.3  Calibrated  HSPF  Model  Parameters  for  the  Athabasca  and  Beaver  River  Basins 


WELL  DRAINED 

TILL 

WELL  DRAINED 

SAND 

WELL  DRAINED 

CLAY  LOAM 

POORLY  DRAINED  TILL 
(LOWLAND 
GLACiOFLUVIAL) 

POORLY  DRAINED  SAND 
(LOWLAND 
GLACIOLACUSTR1NE) 

POORLY  DRAINED 
CLAY  LOAM 
(LOWLAND 
GLACIAL) 

ORGANIC 

IM  PERVIOUS 
(FRACTURED  ROCK 
TREATED  AS 
PERVIOUS) 

Water 

Meaning  

Parameter 

FOREST 

The  fraction  of  the  pervious  Hand  segment  w  hich  is 

none 

0.8 

0.5 

0.8 

0.8 

0.8 

0.8 

0.8 

0.2 

LZSN 

The  lower  zone  nominal  storage 

in 

0.3 

2 

0.05 

0.3 

0.4 

0.9 

13.26 

An  index  to  the  infiltration  capacity  of  the  soil 

in/hr 

0.5 

0.01  73 

0.05 

0.008 

!  o.oi 

0.5 

0.02  (0.04) 

KV  ARY 

parameter  w  hich  affects  the  behavior  of 
groundw  ater  recession  flow  ,  enabling  it  to  be  non¬ 
exponential  in  its  decay  w  ith  time 

1/in 

0.03 

5 

1.18 

o 

O 

o 

2.847 

0.8 

AGWRC 

The  basic  groundwater  recession  rate  if  KV  ARY  is 

zero  and  there  is  no  inflow  to  groundw  ater. 

1/day 

0.993  (0.983) 

0.8 

0.938 

0.87 

0.87 

I  0.87 

0.992 

0.997 

F’ETMAX 

The  air  temperature  below  w  hich  E-T  w  ill  arbitrarily 

be  reduced  below  the  value  obtained  from  the  input 
time  series. 

degF 

40 

40 

40 

40 

40 

j  40 

40 

40 

PETMIN 

The  temperature  below  w  hich  E-T  w  ill  be  zero 

regardless  of  the  value  in  the  input  time  series. 

degF 

35 

35 

35 

35 

35 

!  35 

35 

35 

INFEXP 

Exponent  in  the  infiltration  equation 

none 

2 

2 

2 

2 

2 

1  2 

2 

2 

INFILD 

Ratio  between  the  maximum  and  mean  infiltration 

capacities 

none 

2 

2 

2 

2 

2 

i  2 

2 

2 

DEEPFR 

Fraction  of  groundw  ater  inflow  w  hich  w  ill  enter 

deep  (inactive)  groundwater 

none 

O 

O 

O 

0.1  1 

O 

O 

O 

O 

BASETP 

Fraction  of  remaining  potential  E-T  w  hich  can  be 

satisfied  from  baseflow  (groundwater  outflow  ),  if 
enough  is  available. 

none 

0.005 

0.005 

0.005 

0.3 

0.2 

0.2 

0.005 

0.005 

AGWETP 

Fraction  of  remaining  potential  E-T  w  hich  can  be 
satisfied  from  active  groundwater  storage  if  enough 
is  available. 

none 

0.01 

0.01 

0.01 

0.4 

0.01 

0.01 

0.01 

0.01 

CEPSC 

Interception  storage  capacity. 

in 

see  monthly  table 

see  monthly  table 

see  monthly  table 

see  monthly  table 

0.1 

0.10 

see  monthly  table 

see  monthly  table 

UZSN 

Upper  zone  nominal  storage. 

in 

0.1  (0.2) 

0.05 

0.3 

0.5 

0.5 

0.5 

0.703 

0.6 

NS  UR 

Manning's  n  for  the  overland  flow  plane. 

complex 

0.25 

0.25 

0.25 

0.35 

0.35 

0.35 

0.25 

0.25 

INTFW 

Interflow  inflow  parameter. 

none 

3.3 

4.83 

1 

25 

8 

io 

8.42 

3.3 

IRC 

Interflow  recession  parameter 

1/day 

0.94 

0.798 

0.534 

0.92 

0.925 

0.925 

0.944 

0.2 

LZETH 

Low  er  zone  E-T  parameter. 

none 

see  monthly  table 

see  monthly  table 

see  monthly  table 

see  monthly  table 

0.5 

0.5 

see  monthly  table 

see  monthly  table 

Interception  Monthly  Table 

WELL  DRAINED 

TILL 

WELL  DRAINED 

SAND 

WELL  DRAINED 

CLAY  LOAM 

POORLY  DRAINED  TILL 
(LOWLAND 
GLACIOFLUVIAL) 

POORLY  DRAINED  SAND 
(LOWLAND 

GL_  AC  IOL  ACUSTRINE) 

POORLY  DRAINED 

CLAY  LOAM 
(LOWLAND 
GLACIAL) 

ORGANIC 

IM  PERVIOUS 

(FRACTURE)  ROCK 
TREATED  AS 
PERVIOUS) 

Jan 

0.5 

0.5 

0.5 

1 

M/A 

N/A 

1 

1 

Feb 

0.5 

0.5 

0.5 

1 

N/A 

N/A 

1 

1 

Mar 

0.1 

0.1 

0.1 

1  .2 

N/A 

N/A 

1 .2 

1.2 

Apr 

0.1 

0.1 

0.1 

0.4 

N/A 

N/A 

0.4 

0.4 

May 

0.05 

0.05 

0.05 

0.05 

N/A 

N/A 

0.05 

0.05 

Jun 

0.1 

0.1 

0.1 

0.1 

N/A 

N/A 

0.1 

0.1 

Jul 

0.05 

0.05 

0.05 

0.05 

N/A 

N/A 

0.05 

0.05 

Aug 

0.35 

0.35 

0.35 

0.35 

N/A 

N/A 

0.35 

0.35 

Sep 

0.4 

0.4 

0.4 

0.4 

N/A 

N/A 

0.4 

0.4 

Oct 

0.4 

0.4 

0.4 

0.4 

N/A 

N/A 

0.4 

0.4 

Nov 

0.4 

0.4 

0.4 

0.4 

N/A 

N/A 

0.4 

0.4 

Dec 

0.4 

0.4 

0.4 

0.4 

N/A 

N/A 

0.4 

0.4 

Lower  Zone  Evapo-trans piration  Monthly  Table 

WELL  DRAINED 

TILL 

WELL  DRAINED 

SAND 

WELL  DRAINED 

CLAY  LOAM 

POORLY  DRAINED  TILL 
(LOWLAND 
GLACIOFLUVIAL) 

POORLY  DRAINED  SAND 
(LOWLAND 

GL  AC  IOL  ACUSTRINE) 

POORLY  DRAINED 

CLAY  LOAM 
(LOWLAND 
GLACIAL) 

ORGANIC 

IM  PERVIOUS 
(FRACTURE)  ROCK 
TREATS)  AS 
PERVIOUS) 

Jan 

0.3 

0.3 

0.3 

0.01 

N/A 

N/A 

0.3 

0.3 

Feb 

0.5 

0.5 

0.5 

0.01 

N/A 

N/A 

0.5 

0.5 

Mar 

0.6 

0.6 

0.6 

0.01 

N/A 

N/A 

0.6 

0.6 

Apr 

0.8 

0.8 

0.8 

0.1  ~ * 

N/A 

N/A 

0.8 

os 

M  ay 

0.2 

0.2 

0.2 

0.1 

N/A 

N/A 

0.2 

0. 1 

Jun 

0.2 

0.2 

0.2 

0.1 

N/A 

N/A 

0.2 

0.1 

Jul 

0.2 

0.2 

0.2 

0.1 

N/A 

N/A 

0.2 

0.2 

Aug 

0.4 

0.4 

0.4 

0.1 

N/A 

N/A 

0.4 

0.4 

Se  p 

0.5 

0.5 

0.5 

0.1 

N/A 

N/A 

0.5 

0.5 

Oct 

0.5 

0.5 

0.5 

0.1 

N/A 

N/A 

0.5 

- M - 

Nov 

0.5 

0.5 

0.5 

0.1 

N/A 

N/A 

0.5 

0.2  1 

Dec 

0.6 

0.6 

0.6 

0.01 

N/A 

N/A 

0.6 

02 
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HYDRO-CLIMATE  MODEL  SELECTION  AND  APPLICATION  ON  THE  ATHABASCA  AND  BEAVER  RIVER  BASINS 


Table  4.3  Calibrated  HSPF  Model  Parameters  for  the  Athabasca  and  Beaver  River  Basins  (continued) 


Table  4.3b  HSPF  Model  Parameters  (Water)  -  IMPLND 


GLACIER 

Water 

Meaning 

Units 

Parameter 

NSUR 

Manning's  n  for  the  overland  flow  plane 

none 

1 

RETSC 

The  retention  (interception)  storage  capacity  of  the  surface. 

in 

0 

PETMAX 

The  air  temperature  below  w  hich  E-T  w  ill  arbitrarily  be 
reduced  below  the  value  obtained  from  the  input  time  series. 

degF 

48 

PETMIN 

The  terrperature  below  w  hich  E-T  w  ill  be  zero  regardless  of 

the  value  in  the  input  time  series. 

degF 

40 

RETS 

The  initial  retention  storage. 

in 

0.001 

SURS 

The  initial  surface  (overland  flow )  storage. 

in 

0.001 

Table  4.3c  HSPF  Model  Parameters  (Snow) 


WELL 

DRAINED 

TILL 

WELL 

DRAINED 

SAND 

WELL 

DRAINED 

CLAY  LOAM 

POORLY  DRAINED  TILL 
(LOWLAND 
GLACIOFLUVIAL) 

POORLY  DRAINED  SAND 
(LOWLAND 
GLACIOLACUSTRINE) 

POORLY  DRAINED 

CLAY  LOAM 
(LOWLAND 
GLACIAL) 

ORGANIC 

IMPERVIOUS 
(TREATED  AS 
PERVIOUS) 

GLACIER 

(IMPERVIOUS) 

Snow 

Description 

Units 

Parameter 

LAT 

Latitude 

Degree 

54.3  (57.5) 

54.3  (57.5) 

54.3  (57.5) 

54.3  (57.5) 

54.3  (57.5) 

54.3  (57.5) 

54.3  (57.5) 

54.3  (57.5) 

54.3 

SHADE 

Fraction  of  the  land  w  hich  is  shaded  from  solar 
radiation  by  trees 

none 

0.8 

0.8 

0.8 

0.8 

0.8 

0.8 

0.8 

0.8 

0.66 

SNOWCF 

Factor  by  w  hich  the  input  precipitation  data  w  ill 

be  multiplied 

none 

1 

1 

’  :'1  . 

1 

1 

1 

1 

1 

1 

COVIND 

Maximum  snow  pack  (w  ater  equivalent)  at  w  hich 
the  entire  land  w  ill  be  covered  w  ith  snow 

none 

10 

5 

5 

5 

5 

5 

5 

3 

8.8 

KMELT 

Constant  degree-day  factor  for  the  temperature 

index  snowmelt  method 

in/day.F 

0 

0 

0 

0 

0 

0 

0 

0 

0 

TBASE 

Reference  temperature  for  the  temperature  index 

method 

degF 

32 

32 

32 

32 

32 

32 

32 

32 

32 

RDCSN 

Density  of  cold,  new  snow  relative  to  w  ater 

none 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

TSNOW 

Air  temperature  below  w  hich  precipitation  w  ill  be 
snow 

degF 

40 

40 

40 

40 

40 

40 

37 

40 

30.2 

SNOEVP 

Parameter  w  hich  adapts  the  snow  evaporation 

(sublimation)  equation  to  field  conditions 

none 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0.2 

0.25 

0.0003 

CCFACT 

Parameter  w  hich  adapts  the  snow 
condensation/convection  melt  equation  to  field 
conditions. 

none 

0.1  (0.2) 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1  (0.2) 

0.677 

MW  A  TER 

Maximum  water  content  of  the  snow  pack,  in 

depth  of  water  per  depth  of  water 

none 

0.8 

0.8 

0.8 

0.8 

0.8 

0.8 

0.8 

0.25 

0.005 

MG  MELT 

Maximum  rate  of  snow  melt  by  ground  heat,  in 
depth  of  w  ater  per  day 

in/day 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

0 

0.02 

0 

PACK-ICE 

Uuantity  of  ice  in  the  pack  (w  ater  eauivalentl 

inch 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1000 
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HYDRO-CLIMATE  MODEL  SELECTION  AND  APPLICATION  ON 
THE  ATHABASCA  AND  BEAVER  RIVER  BASINS 


■  Model  parameters  for  the  Well/Rapidly  Drained  Sand  land  type  were  calibrated  using  the  recorded  stream 
flows  at  Driftwood  River  near  the  Mouth  (Environment  Canada  Hydrometric  Station  07BK007).  The  surficial 
geology  of  the  Driftwood  River  sub-basin  is  approximately  56%  Well/Rapidly  Drained  Till,  3%  Organic,  and 
41%  Well/Rapidly  Drained  Sand  (see  Figure  4.2).  During  the  calibration  process,  the  model  parameters  for 
the  Well/Rapidly  Drained  Till  and  Organic  land  types  were  transferred  from  those  obtained  during  the 
calibration  of  the  Lovett  River  and  Mackay  River  above  Dunkirk  River  sub-basins,  respectively.  The 
precipitation  data  used  were  from  the  Slave  Lake  climate  station. 

■  Model  parameters  for  the  Beaver  River  at  Cold  Lake  Reserve  sub-basin  (Environment  Canada  Hydrometric 
Station  06AD006)  were  initially  assumed  to  be  the  same  as  for  the  Well/Rapidly  Drained  Till  land  type  as 
obtained  during  the  calibration  of  the  Lovett  River  sub-basin.  The  entire  Beaver  River  watershed  was 
assumed  to  be  covered  by  Well/Rapidly  Drained  Till  (see  Figure  4.3).  The  calibration  of  the  Beaver  River 
Basin  was  fine  tuned  using  recorded  flows  on  the  basin  because  the  basin  has  a  hydrologic  response  that 
is  significantly  different  from  most  of  the  sub-basins  in  the  Athabasca  River  Basin  due  to  the  significant 
portions  of  the  basin  being  classified  as  non-contributing  during  storms  with  a  frequency  of  more  than  two 
years.  The  precipitation  data  used  were  from  the  Cold  Lake  climate  station. 

■  Model  parameters  for  the  Poorly  Drained  Till  (Lowland  Glaciofluvial),  Poorly  Drained  Sand  (Lowland 
Glaciolacustrine),  and  Poorly  Drained  Clay  Loam  (Lowland  Glacial)  land  types  were  obtained  from  a  report 
on  Regional  Surface  Water  Hydrology  Study  for  Re-Calibration  of  HSPF  Model  (Golder  2003). 

The  sub-basins  used  for  validation  are  listed  in  Table  4.2  and  shown  in  Figure  4.4.  The  validation  of  the 
calibrated  HSPF  model  was  based  on  a  comparison  of  observed  and  simulated  flows  at  gauging  stations  on  the 
main  stem  of  the  Athabasca  River,  namely,  Athabasca  River  near  Jasper  (Environment  Canada  Hydrometric 
Station  07AA002),  Athabasca  River  at  Hinton  (Environment  Canada  Hydrometric  Station  07AD002),  Athabasca 
River  near  Windfall  (Environment  Canada  Hydrometric  Station  07AE001),  Athabasca  River  at  Athabasca 
(Environment  Canada  Hydrometric  Station  07BE001),  Athabasca  River  below  McMurray  (Environment  Canada 
Hydrometric  Station  07DA001),  and  Clearwater  River  at  Draper  (Environment  Canada  Hydrometric  Station 
07CD001).  Model  parameters  for  the  Clearwater  River  at  Draper  (Environment  Canada  Hydrometric  Station 
07CD001)  were  calibrated  and  validated  separately  as  most  of  the  sub-basin  lies  within  the  province  of 
Saskatchewan  and  the  land  type  information  was  not  available.  The  portion  of  the  Clearwater  River  sub-basin 
with  missing  surficial  geology  data  was  assumed  to  be  covered  by  rapidly  drained  sand  (i.e.,  extending  the 
surficial  geology  data  available  on  the  province  of  Alberta  side  of  the  border),  and  the  calibration/validation  was 
done  using  recorded  stream  flow  data. 

The  accuracy  of  the  model  calibration  and  validation  was  evaluated  by  comparing  the  measured  and  simulated 
flow  parameters  listed  below: 

■  Mean  annual  flow; 

■  Mean  open  water  (March  to  October)  flow; 

■  Mean  monthly  flow  (12  months);  and, 

■  2-,  10-,  and  25-year  peak  flood  flows.  The  2,  10,  and  25-yr  flood  flows  were  estimated  from  the  “best-fit” 
probability  distribution  function  to  the  annual  peak  flow  series.  The  available  probability  distribution 
functions  were  the  Three-Parameter  Lognormal,  General  Extreme  Value,  and  Log  Pearson  III 
distributions.  The  “best  fit”  distribution  was  based  on  the  minimum  sum  of  squares  of  differences 
between  the  data  and  the  theoretical  fits. 
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4. 2.1. 5  Model  Calibration  and  Validation  Statistics 

A  comparison  of  the  statistics  of  the  simulated  flow  series  from  the  calibration  and  validation  of  HSPF  with  the 
statistics  of  the  observed  series  is  provided  in  Appendix  E.  Observed  and  simulated  annual  hydrographs  are 
provided  therein  for  selected  time  periods.  Comparisons  of  simulated  and  observed  mean  monthly  flows  on  the 
six  gauged  sub-basins  used  for  calibration,  each  dominated  by  one  particular  land  type,  are  shown  in  Figure  4.5. 
Figure  4.6  shows  the  mean  monthly  flow  comparisons  for  the  seven  sub-basins  used  for  validation  of  the 
calibrated  HSPF  model. 

The  outcome  of  the  calibration  and  validation  of  the  HSPF  model  was  deemed  to  be  good  to  reasonable  to  poor 
based  on  the  following  criteria: 


■  Good: 


■  Observed  mean  annual  flow  or  mean  open-water  flow  replicated  to  less  than  10%;  and, 

■  Mean  monthly  flows  replicated  to  within  20%,  except  for  winter  (very  low  flow)  months  when  a 
difference  of  less  than  40%  was  deemed  to  be  good. 

■  Reasonable: 

■  Observed  mean  annual  flow  or  mean  open-water  flow  replicated  to  less  than  15%;  and, 

■  Mean  monthly  flows  replicated  to  within  40%,  except  for  winter  (very  low  flow)  months  when  a 
difference  of  less  than  60%  was  deemed  to  be  reasonable. 


■  Poor: 


■  Difference  between  observed  and  simulated  mean  annual  flows  is  more  than  15%;  and, 

■  Difference  between  observed  and  simulated  mean  monthly  flows  is  greater  than  40%,  and  greater  than 
60%  for  winter  (very  low  flow)  months. 
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Note:  [1]  Open  water  season  is  from  April  to  October. 
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Note:  [1]  Open  water  season  is  from  May  to  October. 

[2]  There  is  no  recorded  flow  data  during  winter  season. 
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Note:  [1]  Open  water  season  is  from  April  to  October. 
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The  calibration  on  mean  annual  and/or  mean  open-water  flows  is  generally  good  for  five  calibration  sub-basins 
and  poor  on  the  Lovett  River  sub-basin.  With  respect  to  mean  monthly  flows,  the  calibration  is  generally 
reasonable,  although  there  are  some  significant  differences  between  observed  and  simulated  values  for  some 
winter  months,  which  is  not  unexpected,  and  generally  for  the  Miette  River  sub-basin.  The  reasons  for  some  of 
the  significant  variances  between  observed  and  simulated  statistics  can  be  explained  as  follows.  For  the 
calibration  of  these  sub-basins,  climate  data  are  not  generally  available  within  the  sub-basin  itself,  instead,  the 
data  are  transferred  from  other  locations.  For  example,  the  calibration  on  the  Driftwood  River  sub-basin  is  based 
on  climate  data  recorded  at  Slave  Lake  Station  (more  than  40  km  from  the  center  of  the  sub-basin).  The 
drainage  areas  of  the  calibration  sub-basins  are  relatively  small,  and,  therefore,  are  more  likely  to  be  subject  to 
uncertainties  due  to  the  spatial  variability  in  precipitation  than  larger  basins.  The  timing  and  magnitude  of  actual 
within-basin  precipitation  can  be  different  from  the  recorded  data  at  the  climate  station  used  for  calibration. 
Redistribution  of  snow  on  the  landscape  can  have  a  significant  influence  on  the  rate  and  timing  of  snow  melt  and 
the  soil  water  regime,  and  hence  watershed  yield.  These  processes  are  difficult  to  model,  primarily  because  of  a 
general  lack  of  the  required  data,  except  perhaps  in  research  basins.  Hence,  it  is  difficult  to  get  good  calibration 
for  small  sub-basins  that  are  more  prone  to  be  affected  by  localized  precipitation.  In  addition,  climate  stations 
tend  to  be  located  at  relatively  low  elevations  compared  to  the  runoff-producing  areas  in  the  upper  Athabasca 
River  Basin.  Extrapolating  the  station  data  to  high  elevation  sub-basins  can  be  problematic.  For  example, 
precipitation  data  at  the  Jasper  climate  station,  which  is  located  at  an  elevation  of  about  1,050  metres  above  sea 
level  (masl),  was  used  for  simulation  of  runoff  from  the  Miette  River  sub-basin.  The  summer  and  winter 
precipitation  was  adjusted  using  an  orographic  adjustment  factor.  More  than  70%  of  the  Miette  River  sub-basin 
is  at  an  elevation  greater  than  2,000  masl,  with  a  maximum  elevation  up  to  2,500  masl.  It  is  possible  that 
extrapolation  of  the  orographic  adjustment  factor  to  the  very  high  elevations  in  this  sub-basin  may  have  resulted 
in  some  errors  in  the  simulated  runoff  values.  In  practice,  there  should  be  a  cap  on  the  application  of  the 
orographic  adjustment  with  elevation  since  rainfall  amounts  tend  to  decrease  past  certain  elevations. 

The  calibration  shows  significant  differences  between  observed  and  simulated  values  for  winter  flows  in  the 
small  sub-basins  because  (1)  there  are  significant  uncertainties  in  low  winter  flow  values  and  (2)  small 
differences  in  magnitude  usually  manifest  themselves  as  large  percentage  changes  because  the  differences  are 
divided  by  small  winter  flow  values. 

Other  factors  that  may  have  affected  the  calibration  of  the  model  are  spatial  variability  in  frozen  soil  conditions, 
spatial  variability  in  vegetation  cover,  land  use  changes  over  time,  water  withdrawals  and  returns,  etc. 
Adjustments  for  these  factors  can  be  made  within  the  model  for  specific  studies  on  stand-alone  sub-basins, 
however,  given  that  the  model  is  being  implemented  for  the  entire  Athabasca  River  Basin  with  the  focus  on  the 
natural  flows  in  the  lower  reaches  of  the  Athabasca  River  Basin,  these  adjustments  are  not  considered 
necessary  at  this  stage. 

Further  attempts  at  refining  the  calibration  of  the  model  did  not  result  in  significant  improvements  in  model 
performance.  Also,  given  that  the  performance  of  the  model  at  the  validation  nodes  on  the  main  stem  of  the 
Athabasca  River  is  considered  good,  the  HSPF  model  is  considered  calibrated  with  parameters  as  given  in 
Table  4.3  for  the  Athabasca  and  Beaver  River  basins. 

Figure  4.6  and  the  detailed  results  presented  in  Appendix  E  indicate  that,  in  general,  the  calibrated  model 
reproduced  the  measured  discharges  at  the  validation  nodes  on  the  main  stem  of  the  Athabasca  River  well  to 
reasonably  well,  giving  confidence  in  the  use  of  the  model  for  assessing  the  hydrologic  effects  of  potential  future 
climate  changes.  The  validation  nodes  capture  sub-basins  with  different  land  types  (in  different  percentages 
within  each  sub-basin  as  shown  in  Table  4.2)  and  different  climate  regime.  The  combination  of  a  range  of  land 
types  and  climate  regimes  at  the  validation  nodes  is  likely  a  more  rigorous  test  of  the  performance  of  the 
calibrated  model.  However,  it  may  be  argued  that  at  these  nodes  the  drainage  areas  are  much  larger  than  those 
of  the  small  sub-basins  and  differences  in  responses  between  small  sub-basins  tend  to  be  “masked”  out,  thus 
improving  the  model  performance  at  the  validation  stage.  Notwithstanding  the  foregoing,  it  is  concluded  that  the 
calibrated  HSPF  model  has  been  validated  and  is  appropriate  for  assessing  the  effects  of  climate  change  on  the 
yield  in  the  Lower  Athabasca  Regional  Plan  Area. 
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4.3  Summary 

The  HSPF  model  was  calibrated  and  validated  for  application  to  the  Athabasca  River  Basin  and  Beaver  River 
Basin.  Each  basin  was  sub-divided  on  the  basis  of  drainage  network,  locations  of  gauging  stations,  and  surficial 
geology.  Temperature  and  precipitation  data  for  the  calibration  of  the  model  were  from  climate  stations  closest 
(specific  selection  based  on  availability  of  concurrent  climate  and  flow  data)  to  the  sub-basins.  Seven  sub¬ 
basins  were  selected  for  calibration  of  the  model  and  six  (located  primarily  on  the  main  stem  of  the  Athabasca 
River)  were  selected  for  validation  of  the  calibrated  parameters.  The  calibration  on  mean  annual  and/or  mean 
open-water  flows  is  generally  good  for  most  calibration  sub-basins.  With  respect  to  mean  monthly  flows,  the 
calibration  is  generally  reasonable,  although  there  are  some  significant  differences  between  observed  and 
simulated  values  for  some  winter  months.  The  calibrated  model  reproduced  the  measured  discharges  at  the 
validation  nodes  on  the  main  stem  of  the  Athabasca  River  well  to  reasonably  well,  giving  confidence  in  the  use  of 
the  model  for  assessing  the  hydrologic  effects  of  potential  future  climate  changes.  It  was  concluded  that  the 
calibrated  HSPF  model  has  been  validated  and  is  appropriate  for  assessing  the  effects  of  climate  change  on  the 
yield  in  the  Lower  Athabasca  Regional  Plan  Area. 
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5.0  ASSESSMENT  OF  HYDROLOGIC  EFFECTS  OF  CLIMATE  CHANGE 

5.1  Introduction 

The  effects  of  climate  change  on  water  yield  from  the  basins  of  Alberta  will  affect  water  uses  and  water 
management  in  those  basins.  There  is  a  need  to  assess  the  potential  effects  so  that  watershed  planners  can 
adapt  their  plans  to  take  advantage  of  positive  effects  and  implement  mitigation  measures  to  minimize  the 
negative  effects.  Alberta  is  embarking  on  a  regional  planning  process  and  the  effect  of  climate  change  on 
watershed  responses  is  a  key  consideration  in  that  process. 

5.2  Baseline  Climate  Conditions  and  Future  Climate  Scenarios 

To  address  the  potential  effects  of  climate  change  on  water  yield  in  the  Lower  Athabasca  Regional  Plan  Area 
(LARP),  which  includes  the  downstream  portion  of  the  Athabasca  River  Basin  and  the  Beaver  River  Basin,  the 
HSPF  model,  calibrated  and  validated  for  these  two  basins,  has  been  used  to  simulate  the  hydrologic  effects  of 
forecasted  future  climate  scenarios.  An  analysis  of  the  effects  of  climate  change  depends  not  only  on  future 
conditions  but  also  on  the  baseline  climate  to  which  the  predictions  are  compared.  The  IPCC  recommends  that 
1961  to  1990  be  adopted  as  the  climatological  baseline  period  in  impact  assessments. 

As  a  first  step  in  addressing  the  potential  effects  of  climate  change  on  key  ecosystem  variables  in  Alberta, 
Alberta  Environment  has  developed  a  comprehensive  database  (Alberta  Climate  Model)  of  thirteen  climate 
variables  (including  annual  mean  temperature,  annual  maximum  temperature,  annual  minimum  temperature, 
annual  precipitation,  degree-days  above  a  threshold,  annual  moisture  index)  from  the  available  Environment 
Canada’s  climate  records  from  1961  to  1990  (Alberta  Environment  2005).  The  Alberta  Climate  Model  combines 
long-term  records  from  weather  stations  with  physiographic  descriptors  of  station  locations  to  allow  the 
interpolation  and  extrapolation  of  values  for  climatic  variables  along  gradients.  Using  mathematical  techniques, 
the  pattern  of  fixed-point  climate  stations  has  been  transformed  into  continuous  surfaces  that  allow  the 
estimation  of  climatic  variables  at  any  point  in  Alberta.  The  Alberta  Climate  Model  is  being  used  in  a  variety  of 
research  projects  (including  reclamation  and  terrestrial  vegetation).  The  Alberta  Climate  Model  facilitates 
estimation  of  values  for  climatic  variables  in  areas  poorly  represented  by  climatic  stations  as  well  as  provides  a 
baseline  from  which  to  estimate  impacts  of  changes  suggested  for  future  climates. 

The  Alberta  Climate  Model  is  not  a  tool  for  predicting  future  climates,  but  it  contains  predictions  for  future 
climates.  Barrow  and  Yu  (2005)  combined  the  climate  information  derived  from  the  Alberta  Climate  Model  with 
changes  in  climate  variables  predicted  by  a  selected  number  of  GCMs-scenarios  to  provide  “actual  future” 
climate  scenarios  for  the  2020s,  2050s  and  2080s.  The  sub-set  of  climate  change  scenarios  was  selected  on 
the  basis  of  changes  in  summer  mean  temperature  and  precipitation  for  the  2050s.  The  five  selected  scenarios 
represented  conditions  that  were  cooler  and  wetter  (NCARPCM  A1B),  cooler  and  drier  (CGCM2  B2(3)),  warmer 
and  wetter  (HadCMS  A2(a))  and  warmer  and  drier  (CCSRNIES  A1  FI)  than  median  conditions  (HadCM3  B2(b)). 
These  five  scenarios  bound  the  range  of  possible  future  climates  from  a  number  of  GCMs  with  respect  to 
temperature  and  precipitation. 

Because  of  the  differences  in  scale  between  the  GCM  grids  (2.5°  to  5.625°)  and  the  resolution  (1  km)  of  the 
Alberta  Climate  Model,  Barrow  and  Yu  (2005)  interpolated  the  climate  change  scenarios  to 
0.5°  latitude/longitude  resolution  using  a  bilinear  two-dimensional  interpolation  routine.  The  interpolated  climate 
change  scenario  was  then  applied  to  the  observed  baseline  climatology.  Each  0.5°  resolution  scenario  grid  box 
contains  60  grid  boxes  for  the  baseline  climatology. 

AENV  provided  Golder  with  the  Alberta  Climate  Model  data  and  average  monthly  temperature  and  precipitation 
predicted  by  five  GCMs  for  two  future  periods:  2010  to  2039  (referred  to  as  2020s)  and  2040  to  2069  (referred  to 
as  2050s).  An  assessment  of  the  effects  of  the  2080s  scenarios  was  not  within  the  scope  of  this  study.  The 
Alberta  Climate  Model  provides  the  average  monthly  temperature  and  precipitation  and  represents  baseline 
conditions  for  all  areas  of  Alberta  from  1961  to  1990.  Average  (1961-1990)  monthly  temperature  and 
precipitation  values  are  made  available  by  AENV  for  grids  of  size  of  approximately  1  km  by  1  km  (0.0083333°) 
covering  the  entire  province. 
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The  five  climate  scenarios  (2020s  and  2050s)  provided  by  AENV  were: 

m  CCSRNIES_A1  FI  (warmer  and  drier  than  median  conditions); 
h  CGCM2_B23  (cooler  and  drier  than  median  conditions); 
m  HADCM3_A2A  (warmer  and  wetter  than  median  conditions); 

■  HADCM3J32B  (median  conditions);  and, 

m  NCARPPCM_A1  B  (cooler  and  wetter  than  median  conditions). 

The  average  (2010-2039  and  2040-2069)  monthly  temperature  and  precipitation  values  are  made  available  for 
the  same  grids  as  the  baseline  climate  for  the  entire  province. 

5.3  Assessment  of  Climate  Change  Effects  -  Athabasca  and  Beaver 
River  Basins 

5.3.1  Baseline  Climate  for  Model  Simulations 

The  Alberta  Climate  Model  provides  only  the  average  (1961-1990)  monthly  temperature  and  precipitation  values. 
The  HSPF  model  is  a  continuous  simulation  model  that  requires  temperature  and  precipitation  inputs  as  daily 
values.  The  baseline  input  daily  series  for  HSPF  were  obtained  from  the  records  at  climate  stations  within  the 
Athabasca  and  Beaver  River  basins.  Using  the  same  data  from  which  the  Alberta  Climate  Model  was  derived, 
daily  series  of  temperature  and  precipitation  were  compiled  for  six  index  climate  stations  (see  map  in  Appendix 
D)  in  the  Athabasca  River  Basin.  The  criteria  used  to  select  the  index  stations  were:  range  of  climate  variables 
recorded  by  the  station  (e.g.,  precipitation,  temperature,  wind  speed,  etc),  length  of  recorded  data,  preferable 
from  1960  to  the  present  and  at  least  covering  the  baseline  period  of  1961  to  1990,  reliable/good  quality  data 
with  low  number  of  missing  data,  and  spatial  distribution  to  cover  most  of  the  Hydrologic  Regions  in 
Athabasca/Beaver  river  basins.  The  index  climate  stations  are: 

■  Jasper  Warden,  covering  sub-basins  in  Hydrologic  Regions  3  and  4; 

■  Edson,  covering  sub-basins  in  Hydrologic  Regions  5  and  10  South; 

■  Campsie,  covering  sub-basins  in  Hydrologic  Region  8  South; 

■  Slave  Lake,  covering  sub-basins  in  Hydrologic  Regions  10  North  and  8  North; 

■  Lac  La  Biche,  covering  sub-basins  in  Hydrologic  Region  2C;  and, 

■  Fort  McMurray,  covering  sub-basins  in  Hydrologic  Region  9A. 

These  stations  provide  good  coverage  of  the  Athabasca  River  Basin.  The  Cold  Lake  climate  station  was  used 
as  the  index  station  for  the  Beaver  River  Basin. 

Missing  data  at  the  index  stations  were  filled  in  by  first  developing  a  relationship  between  seasonal  precipitation 
and  elevation  using  data  at  all  climate  stations  in  a  Hydrologic  Region,  and  then  using  the  relationship  to  transfer 
climate  data  from  nearby  climate  stations  adjusted  for  elevation  to  the  index  station  within  the  same  Hydrologic 
Region. 

5.3.2  Baseline  Climate  for  Estimating  Changes  in  Future  Temperature  and 

Precipitation 

The  baseline  (1961-1990)  average  monthly  temperature  and  precipitation  values  for  each  of  the  sub-basins 
included  in  the  HSPF  model  were  estimated  as  the  average  of  the  grid  cell  values  within  the  sub-basins.  The 
average  data  from  each  sub-basin  was  then  aggregated  to  the  corresponding  “Hydrologic  Region”  and  applied  to 
the  Index  Climate  Station.  Elevation  correction  factors  are  then  applied  as  necessary. 

The  results  are  provided  in  Appendix  F. 
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5.3.3  Estimating  Changes  in  Future  Temperature  and  Precipitation 

The  2020s  and  2050s  average  monthly  temperature  and  precipitation  values  for  each  of  the  sub-basins  included 
in  the  HSPF  model  were  estimated  as  the  average  of  the  grid  cell  values  within  the  sub-basins.  The  average 
data  from  each  sub-basin  was  then  aggregated  to  the  corresponding  “Hydrologic  Region”  and  applied  to  the 
Index  Climate  Station.  Elevation  correction  factors  are  then  applied  as  necessary.  The  differences  between  the 
baseline  averages  and  the  2020s  averages,  and  the  differences  between  the  baseline  averages  and  the  2050s 
averages  are  shown  in  Appendices  G  and  H,  respectively,  for  each  sub-basin  in  the  Athabasca  River  and 
Beaver  River  Basins.  The  differences  in  precipitation  are  shown  as  the  ratio  of  the  future  value  to  the  baseline 
value,  and  the  differences  in  temperature  are  shown  as  the  future  value  minus  the  baseline  value. 

Table  5.1  summarizes  the  changes  in  precipitation  and  temperature  as  an  average  for  the  entire  Athabasca 
River  Basin  for  each  climate  change  scenario  for  the  2020s  and  2050s.  In  a  very  few  instances,  for  example, 
change  in  precipitation  for  January  using  the  CCSRNIES-A1F1  model/scenario,  Table  5.1  shows  the  value  as 
“missing”.  The  software  for  processing  the  Climate  Model  data  could  not  read  the  block  of  data  for  the  particular 
month  and  climate  model/scenario.  It  was  not  immediately  apparent  after  some  investigation  where  the  source 
of  the  problem  was  within  the  block  of  data  for  the  particular  month  (1680  rows  x  1800  columns  for  each  variable 
for  each  month).  Because  of  their  relative  rare  occurrences,  the  problem  data  block  was  set  aside.  For  further 
modelling  work,  the  “missing”  value  was  replaced  as  the  average  of  the  available  values  for  the  two  adjacent 
months.  Figures  5.1  and  5.2  show  the  changes  in  mean  annual  precipitation  and  mean  annual  temperature, 
respectively,  for  the  2020s  and  2050s  compared  to  the  baseline  values. 
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HYDRO-CLIMATE  MODEL  SELECTION  AND  APPLICATION  ON 
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Figure  5. 1  Percent  Change  in  Precipitation  for  2020s  and  2050s  -  Athabasca  River  Basin 
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Figure  5.2  Change  in  Mean  Annual  Temperature  for  2020s  and  2050s  -  Athabasca  River  Basin 
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Table  5.1  shows  that  the  average  change  in  mean  annual  precipitation  varies  from  -1%  to  +7%  for  the  2020s 
scenarios  and  from  +2%  to  +13%  for  the  2050s  scenarios.  The  range  of  the  change  in  precipitation  is  much 
wider  on  a  monthly  basis:  -10%  to  +23%  for  the  2020s  scenarios  and  -15%  to  +46%  for  the  2050s  scenarios. 
The  general  trend  appears  to  be  increased  precipitation  from  the  baseline  period  to  the  2050s,  but  with  greater 
variability  in  the  monthly  changes.  Table  5.1  shows  that  the  average  increase  in  mean  annual  temperature 
varies  from  0.63°C  to  1.52°C  for  the  2020s  scenarios  and  from  1.77°C  to  4.35°C  for  the  2050s  scenarios.  The 
range  of  the  change  in  temperature  depends  on  the  season  under  consideration:  increases  are  generally  greater 
for  the  fall  and  winter  months  compared  to  the  summer  months.  The  patterns  of  basin-wide  annual  changes  due 
to  the  five  scenarios  do  not  appear  to  correspond  precisely  with  the  Barrow  and  Yu  (2005)  classifications,  that  is, 
CCSRNIES  A1F1:  warmer  and  drier;  CGCM2  B2(3):  cooler  and  drier;  HadCM3  A2A:  warmer  and  wetter; 
HadCM3  B2B:  median  conditions;  and  NCARPCM  A1B:  cooler  and  wetter.  It  is  noted  that  the  sub-set  of  climate 
change  scenarios  was  selected  by  Barrow  and  Yu  (2005)  on  the  basis  of  changes  in  summer  mean  temperature 
and  precipitation  for  the  2050s.  Nevertheless,  the  five  selected  scenarios  represent  a  range  of  precipitation  and 
temperature  conditions  that  are  useful  for  planning  purposes. 

Table  5.2  summarizes  the  basin  average  changes  in  precipitation  and  temperature  for  the  Beaver  River  Basin  for 
each  climate  change  scenario  for  the  2020s  and  2050s.  The  basin  average  changes  in  mean  annual 
precipitation  and  temperature  in  the  Beaver  River  Basin  are  generally  very  similar  to  those  for  the  Athabasca 
River  Basin  (Table  5.1),  except  that  temperature  increases  appear  to  be  higher  by  between  0.1  °C  and  0.5°C  in 
the  Beaver  River  Basin. 
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5.3.4  Future  Daily  Climate  Scenarios 

Daily  30-year  series  of  temperature  and  precipitation  for  2010-2039  and  2040-2069  were  then  generated  by 
adjusting  the  baseline  climate  data  generated  as  described  in  Section  5.3.1  by  the  differences  estimated  as 
described  in  Section  5.3.3. 

5.3.5  Hydrologic  Effects  of  Climate  Change  on  Flows  in  the  Athabasca  and  Beaver 
River  Basins 

The  HSPF  model  calibrated  for  the  Athabasca  River  Basin  and  Beaver  River  Basin  was  run  with  the  baseline 
climate  data  (Section  5.3.1)  and  the  adjusted  future  climate  data  (Section  5.3.4).  A  comparison  of  the  hydrologic 
statistics  (annual  and  monthly  water  yield  particularly)  at  selected  locations  on  the  main  stem  of  the  Athabasca 
River  and  on  the  Beaver  River  at  the  Cold  Lake  Reserve  is  provided  in  Appendix  I  for  the  2020s  climate 
scenarios  and  in  Appendix  J  for  the  2050s  climate  scenarios.  Summary  statistics  (mean  annual  flow,  mean 
seasonal  flows,  2-yr,  10-yr,  and  25-yr  flood  flows,  and  mean  monthly  flows)  for  the  2020s  and  2050s  climate 
scenarios  are  presented  in  Tables  5.3  and  5.4,  respectively.  The  forecasted  changes  in  mean  monthly  flows  for 
the  2020s  and  2050s  scenarios  are  presented  graphically  in  Figure  5.2  and  Figure  5.3,  respectively.  The 
changes  in  Figure  5.2  and  Figure  5.3  are  shown  in  terms  of  both  percent  changes  and  absolute  flows. 

Tables  5.3  and  5.4  indicate  that  the  hydrologic  effects  of  forecasted  changes  in  climate  vary  both  spatially  across 
the  Athabasca  River  Basin  and  temporally  during  the  year.  General  conclusions  from  the  summary  statistics  are 
as  follows: 

■  All  the  five  climate  scenarios  for  the  2020s  tend  to  result  in  lower  mean  annual  flows,  with  the  decrease  in 
mean  annual  flow  becoming  more  severe  downstream  along  the  Athabasca  River.  Among  the  five 
scenarios,  the  CCSRNIES-A1F1  scenario  predicts  the  largest  decrease,  ranging  from  -8%  near  Jasper  to  - 
16%  near  Fort  McMurray.  For  the  Beaver  River  Basin,  the  CCSRNIES-A1F1  scenario  predicts  a  decrease 
in  mean  annual  flow  of  almost  22%.  The  hydrologic  response  to  the  effects  forecasted  by  the  CCSRNIES- 
A1F1  scenario  appears  to  be  consistent  with  its  “warmer  and  drier  than  median  conditions”  classification  by 
Barlow  and  Yu  (2005).  The  FIADCM3-A2A  scenario  predicts  the  smallest  decrease,  ranging  from  -2%  near 
Jasper  to  -5%  near  Athabasca  and  -2%  near  Fort  McMurray.  For  the  Beaver  River  Basin,  the  HADCM3- 
A2A  scenario  predicts  no  change  in  mean  annual  flow.  The  hydrologic  response  to  the  effects  forecasted 
by  the  HADCM3-A2A  scenario  appears  to  be  somewhat  consistent  with  its  “warmer  and  wetter  than  median 
conditions”  classification  by  Barlow  and  Yu  (2005). 

■  Changes  in  mean  monthly  flows  for  the  2020s  scenarios  tend  to  be  positive  (increase  relative  to  baseline 
values)  during  the  spring  months  and  significantly  negative  (decrease)  during  the  summer  months.  The 
seasonal  differences  in  flows  are  consistent  with  the  relative  effects  of  increased  precipitation  (more  spring 
runoff)  and  increased  temperatures  during  the  summer  months  (higher  summer  evapotranspiration).  The 
differences  in  monthly  hydrologic  responses  between  the  2020s  climate  scenarios  tend  to  be  more  distinct 
in  the  upper  of  the  Athabasca  River  Basin,  and  less  so  near  Fort  McMurray.  It  appears  that  differences  in 
timing  of  runoff  from  the  headwaters  to  the  foothills  sub-basins  tend  to  dampen  the  effects  on  climate 
change  in  the  lower  reaches  of  the  Athabasca  River. 

■  The  2050s  climate  scenarios  result  in  changes  in  mean  annual  flow  that  are  generally  more  severe 
compared  to  the  2020s  scenarios,  which  is  consistent  with  the  changes  in  temperature  and  precipitation 
shown  in  Table  5.1.  Among  the  five  scenarios,  the  CCSRNIES-A1F1  2050s  scenario  predicts  the  greatest 
decrease  in  mean  annual  flows,  ranging  from  -16%  near  Jasper  to  -21%  near  Fort  McMurray.  However,  for 
the  Beaver  River  Basin,  the  CCSRNIES-A1F1  2050s  scenario  predicts  a  decrease  in  mean  annual  flow  of 
9%  compared  to  a  decrease  of  22%  for  the  2020s  scenario.  The  2050s  CCSRNIES-A1  FI  scenario 
predicts  a  higher  increase  in  precipitation  in  the  Beaver  River  Basin  compared  to  the  Athabasca  River 
Basin.  The  HADCM3-B2B,  among  the  five  selected  scenarios,  tends  to  represent  the  median  changes  in 
hydrologic  responses,  which  is  consistent  with  its  classification  by  Barlow  and  Yu  (2005). 
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Table  5.3  Hydrologic  Effects  of  Forecasted  Climate  Change  on  Athabasca  River  and  Beaver  River  Flows  -  2020s  Scenarios 

I - 1  Athabasca  River  near  Jasper  WSC  07AA002  - - 


¥¥ 


g 

I" 


?¥ 


11 


¥?¥ 


I  s  s 


lnm 

i 

1 


¥¥ 


II 

ill 

ill 


si  I 

I  111 


¥¥?¥?88¥ 


S2S 


88¥«8?8«8Sf¥ 


??¥¥ 


! 


r 


■ 


«S 


¥¥ 


II 


ll 

III 

ik 


¥¥¥¥«8 


¥¥¥£¥*¥?¥¥¥¥ 


¥¥««8*¥lS*¥¥ 


¥¥«?¥¥??«?¥¥ 


’Z'smiZznm 


¥88* 


IHilM 


!! 


!i 


ill 


¥¥¥ 


m 


g? 


Ill 


*¥ 


III 

?» 


¥¥¥¥¥¥l'¥8«8 


m 


¥?¥¥S 


Ml 


**¥?¥?¥¥ 


¥¥?« 


8¥??S¥¥ 


*«« 


¥S| 


**? 


ssgg  g 


i 


88 


111 


ill 


88 


88 


¥  I  ¥  ¥ 


5  *  * 


11 


IS 


*¥*88«8fS*¥ 


8f¥?88«¥¥¥¥¥ 


*¥¥¥ 


gggggg 


ggggfggffggg 


^»£g|§§§5Ssg 


gggfg 


gg 


i! 


I 


H 


si 


8«* 


8«! 


|¥¥l¥g 


I! 


n 


n 


m 


in 

I11 


¥¥88¥s||8|8¥ 


f|¥8¥8?|Slp; 


¥¥**88«¥f|f 


C^CNCNJ^CNtCNJCOCOCNJCNCVJCN 


|888fff 


JSi 
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Figure  5.3  Forecasted  Effects  of  Climate  Change  on  Simulated  Mean  Monthly  Flows  -  2020s  Scenarios 
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Figure  5.3  Forecasted  Effects  of  Climate  Change  on  Simulated  Mean  Monthly  Flows  -  2020s  Scenarios 
(continued) 
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Figure  5.3  Forecasted  Effects  of  Climate  Change  on  Simulated  Mean  Monthly  Flows  -  2020s  Scenarios 
(continued) 
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Figure  5.3  Forecasted  Effects  of  Climate  Change  on  Simulated  Mean  Monthly  Flows  -  2020s  Scenarios 
(continued) 
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Figure  5.3  Forecasted  Effects  of  Climate  Change  on  Simulated  Mean  Monthly  Flows  -  2020s  Scenarios 
(continued) 
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Figure  5.4  Forecasted  Effects  of  Climate  Change  on  Simulated  Mean  Monthly  Flows  -  2050s  Scenarios 
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Figure  5.4  Forecasted  Effects  of  Climate  Change  on  Simulated  Mean  Monthly  Flows  -  2050s  Scenarios 
(continued) 
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Figure  5.4  Forecasted  Effects  of  Climate  Change  on  Simulated  Mean  Monthly  Flows  -  2050s  Scenarios 
(continued) 
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Figure  5.4  Forecasted  Effects  of  Climate  Change  on  Simulated  Mean  Monthly  Flows  -  2050s  Scenarios 
(continued) 
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Figure  5.4  Forecasted  Effects  of  Climate  Change  on  Simulated  Mean  Monthly  Flows  -  2050s  Scenarios 
(continued) 
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■  Changes  in  mean  monthly  flows  for  the  2050s  scenarios  tend  to  mimic  the  2020s  results. 

■  The  range  of  effects  of  future  climate  regimes  on  flows  in  the  Lower  Athabasca  Regional  Plan  Area  (LARP) 
is  summarized  in  Table  5.5.  Three  locations  are  used  for  the  summary:  (1)  Athabasca  River  at  Athabasca 
(07BE001)  to  represent  inflow  to  the  LARP  area,  (2)  Athabasca  River  below  McMurray  (07DA001)  to 
represent  flows  in  the  oil  sands  region  where  the  amount  of  water  withdrawal  is  a  key  consideration  in  the 
planning  process,  and  (3)  Beaver  River  at  Cold  Lake  Reserve  sub-basin  (06AD006)  to  represent  the 
southern  portion  of  the  LARP  area.  The  statistics  selected  are  (1)  mean  annual  flow  to  represent  annual 
yield,  (2)  mean  August  flow  to  represent  summer  low  flows,  (3)  mean  February  flow  to  represent  winter  low 
flows,  (4)  mean  June  flow  for  Athabasca  River  and  mean  April  flow  for  Beaver  River  to  represent  spring 
flood  flows,  and  (5)  10-year  flood  flow  to  represent  annual  floods.  For  each  of  the  five  statistics,  Table  5.5 
shows  the  scenario  resulting  in  the  maximum  decrease,  maximum  increase/least  decrease  and  average 
change. 

The  summary  results  indicate  that  water  yield  from  the  two  basins  will  generally  decrease,  with  the  effects 
(in  percentage  terms)  more  significant  during  the  month  of  August.  The  effects  on  low  February  flows, 
although  reduced,  tend  to  be  less  (in  percentage  terms)  compared  to  changes  in  August  flows.  Flood  flows 
tend  to  be  higher  under  the  future  climate  scenarios.  Table  5.5  summarizes  the  effects  of  the  simulated 
climate  scenarios  on  the  10-year  7-day  and  30-day  winter  low  flows  on  the  Athabasca  River  at  Athabasca, 
Athabasca  River  below  McMurray  and  on  the  Beaver  River  at  Cold  Lake.  Tables  5.6,  5.7  and  5.8  provide 
the  winter  low  flow  statistics  at  these  three  locations,  respectively.  The  effects  appear  to  be  very  similar  to 
those  on  the  mean  February  flows  for  the  Athabasca  River  below  Fort  McMurray.  On  the  Beaver  River,  the 
effects  of  the  climate  scenarios  tend  to  be  more  severe  on  the  7-day  and  30-day  lows  than  on  the  mean 
February  flows.  As  expected,  on  a  large  basin,  such  as  the  Athabasca  River  at  Fort  McMurray,  the 
variability  in  winter  low  flows  tend  to  be  much  less  than  on  the  relatively  smaller  basins  such  as  the  Beaver 
River  Basin.  While  the  results  on  the  flood  and  low  flows  are  useful,  it  should  be  noted  that  they  are  based 
on  the  simulation  of  climate  change  characterized  by  predicted  changes  in  mean  monthly  temperature  and 
precipitation.  The  climate  changes  do  not  reflect  changes  that  may  occur  for  shorter  durations  down  to  the 
daily  scale  (scale  at  which  extreme  events  occur),  which  may  influence  the  future  variability  in  the  low  and 
flood  flow  series.  Section  5.4  discusses  the  effect  of  climate  change  on  variability. 
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Table  5.5  Range  of  Potential  Future  Climate  Effects  on  Flows  in  the  Lower  Athabasca  Region 


Flow  Statistic 

Athabasca  River  at  Athabasca 

Athabasca  River  at  McMurray 

Beaver  River  at  Cold  Lake 

Range  of  Change  in  Flow  (%) 

Range  of  Change  in  Flow  (%) 

Range  of  Change  in  Flow  (%) 

2020s 

Low 

High 

Median 

Low 

High 

Median 

Low 

High 

Median 

Mean  Annual  Flow 

-16% 

-5% 

-9% 

-15% 

-2% 

-6% 

-22% 

0% 

-7% 

Scenario 

Cl 

HI 

N 

Cl 

H2 

N 

Cl 

HI 

N  j 

Mean  August  Flow 

-19% 

-11% 

-16% 

-21% 

-6% 

-14% 

-30% 

-2% 

-22% 

Scenario 

Cl 

HI 

C2 

Cl 

HI 

H2 

Cl 

HI 

H2 

Mean  February  Flow 

-15% 

-4% 

-7% 

-12% 

-1% 

-4% 

-26% 

-3% 

-10% 

Scenario 

Cl 

N 

HI 

Cl 

N 

C2 

Cl 

HI 

H2 

Mean  June  or  April  Flow 

-17% 

7% 

-4% 

-17% 

7% 

-5% 

-4% 

11% 

3% 

Scenario 

Cl 

H2 

C2 

Cl 

H2 

N 

C2 

HI 

N 

10-year  Flood  Flow 

-17% 

2% 

-7% 

-18% 

4% 

-3% 

-22% 

17% 

10% 

Scenario 

N 

C2 

H2 

Cl 

HI 

C2 

Cl 

HI 

H2 

10-year  7-day  Low  Flow 

-16% 

-1% 

-4% 

-12% 

2% 

-3% 

-44% 

-7% 

-20% 

Scenario 

Cl 

H2 

HI 

Cl 

H2 

N 

Cl 

C2 

N  ' 

10-year  30-day  Low  Flow 

-16% 

-1% 

-5% 

-12% 

1% 

-3% 

-44% 

-7% 

-20% 

Scenario 

Cl 

H2 

HI 

Cl 

H2 

N 

Cl 

C2 

N 

2050s 

Low 

High 

Median 

Low 

High 

Median 

Low 

High 

Median 

Mean  Annual  Flow 

-21% 

-4% 

-12% 

-21% 

-4% 

-8% 

-19% 

-2% 

-10% 

Scenario 

Cl 

N 

HI 

Cl 

N 

H2 

C2 

HI 

H2 

Mean  August  Flow 

-34% 

-10% 

-24% 

-31% 

-12% 

-21% 

-27% 

-15% 

-24% 

Scenario 

Cl 

N 

C2 

Cl 

N 

HI 

H2 

HI 

N 

Mean  February  Flow 

-17% 

0% 

-8% 

-13% 

1% 

-8% 

-19% 

-7% 

-12% 

Scenario 

Cl 

N 

H2 

Cl 

N 

HI 

C2 

HI 

H2 

Mean  June  or  April  Flow 

-12% 

2% 

-5% 

-15% 

4% 

-3% 

-5% 

18% 

5% 

Scenario 

Cl 

H2 

C2 

Cl 

HI 

N 

C2 

HI 

Cl 

10-year  Flood  Flow 

-11% 

-3% 

-7% 

-15% 

-2% 

-4% 

-11% 

24% 

-3% 

Scenario 

Cl 

N 

C2 

Cl 

HI 

H2 

N 

HI 

Cl 

10-year  7-day  Low  Flow 

-21% 

-2% 

-11% 

-15% 

-2% 

-7% 

-39% 

-20% 

-34% 

Scenario 

Cl 

N 

HI 

Cl 

N 

HI 

C2 

Cl 

H2 

10-year  30-day  Low  Flow 

-22% 

-3% 

-11% 

-16% 

-2% 

-8% 

-39% 

-21% 

-34% 

Scenario 

Cl 

N 

HI 

Cl 

N 

HI 

C2 

Cl 

H2 

Labels  for  Climate  Scenarios 

CCSRNIES-A1F1  Cl 

CGCM-B23  C2 

HADCM3-A2A  |HT~" 

HADCM3-B2B  H2 

NCARPCM-A1B  N 
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Table  5.6  Effects  of  Climate  Scenarios  on  Winter  Low  Flows  in  the  Athabasca  River  at  Athabasca 


Statistic 

Athabasca  River  at  Atl 

habasca  WSC  07BE001 

Baseline  Flows  (1961- 

1990) 

CCSRNIES-A1  FI -2020s 

CGCM-B23-202QS 

HADCM3-A2A-2020S 

HADCM3-B2B-2020S 

NCARPCM-A1  B-2020s 

7Q  10-Low  Flow  (m3/s) 

60.2 

50.3 

54. 

8 

57.6 

59 

.8 

57.2 

30Q  10-Low  Flow  (m3/s) 

63.1 

52.8 

57.2 

|  60.1 

62 

.3 

60.0 

Statistic 

Baseline  Flows  (1961- 

1990) 

CCSRNIES-A1  FI -2050s 

CGCM-B23-2Q50S 

HADCM3-A2A-2050S 

HADCM3-B2B-2050S 

NCARPCM-A1  B-2050s 

7Q  10-Low  Flow  (m3/s) 

60.2 

47.8 

49.8 

53.6 

56 

.7 

58.9 

30Q  10-Low  Flow  (m3/s) 

63.1 

49.3 

51.7 

|  55.6 

59 

.1 

61.4 

Statistic 

Athabasca  River  at  Athabasca  WSC  07BE001 

Baseline  Flows  (1961- 
1990) 

CCSRNIES-A1  FI -2020s 

CGCM-B23-202GS 

HADCM3-A2A-2020S 

HADCM3-B2B-2020S 

NCARPCM-A1  B-2020s 

Change  (%) 

Change  (%) 

Change  (%) 

Change  (%) 

Change  (%) 

7Q1 0-Low  Flow  (m3/s) 

60.2 

-16% 

-8.8% 

-4.2% 

-0.6% 

-5.0% 

30Q1 0-Low  Flow  (m3/s) 

63.1 

-16% 

-9.4% 

-4.7% 

-1.2% 

!  -4.9% 

Statistic 

Baseline  Flows  (1961  - 
1990) 

CCSRNIES-A1  FI -2050s 

CGCM-B23-205GS 

HADCM3-A2A-2050S 

HADCM3-B2B-2050S 

NCARPCM-A1  B-2050s 

Change  (%) 

Change  (%) 

Change  (%) 

Change  (%) 

Change  (%) 

7Q  10-Low  Flow  (m3/s) 

60.2 

-21% 

-17% 

-10.8% 

-5.7% 

-2.1% 

30Q  10-Low  Flow  (m3/s) 

63.1 

-22% 

-18% 

-11.9% 

-6.4% 

-2.7% 

NOTES: 

7Q10  -  10-year  7-day  winter  low  flow  (November/15  ~  April/14) 
30Q10  =  10-year  30-day  winter  low  flow  (November/15  ~  Aprii/14) 
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Table  5.7  Effects  of  Climate  Scenarios  on  Winter  Low  Flows  in  the  Athabasca  River  at  Fort  McMurray 


Athabasca  River  below 

McMurray  WSC  07DA001 

Statistic 

Baseline  Flows  (1961- 

1990) 

CCSRNIES-A1  FI -2020s 

CGCM-B23-2020S 

HADCM3-A2A-2020S 

HADCM3-B2B-2020S 

NCARPCM-A1  B-2020s 

7Q1 0-Low  Flow  (m3/s) 

107 

94.0 

101 

105 

109 

104 

30Q  10-Low  Flow  (m3/s) 

110 

96.7 

104 

108 

112 

107 

Statistic 

Baseline  Flows  (1961- 

1990) 

CCSRNIES-A1  FI -2050s 

CGCM-B23-2O50S 

HADCM3-A2A-2050S 

HADCM3-B2B-2050S 

NCARPCM-A1 B-2050S 

7Q1 0-Low  Flow  (m3/s) 

107 

91.0 

92.8 

!  99.5 

104 

105 

30Q  10-Low  Flow  (m3/s) 

110 

92.8  j 

94.9 

102 

107 

108 

Athabasca  River  below  McMurray  WSC  07DA001 

Baseline  Flows  (1961- 

CCSRNIES-A1  FI -2020s 

CGCM-B23-2020S 

HADCM3-A2A-2020S 

HADCM3-B2B-2020S 

NCARPCM-A1 B-2020S 

Statistic 

1990) 

Change  (%) 

Change  (%) 

Change  (%) 

Change  (%) 

Change  (%) 

7Q1 0-Low  Flow  (m3/s) 

107 

-12% 

-5.4% 

-1.8% 

1.8}% 

-2.8% 

30Q1 0-Low  Flow  (m3/s) 

110 

-12% 

-6.0% 

-2.3% 

1.2% 

-2.7% 

Baseline  Flows  (1961- 

CCSRNIES-A1  FI -2050s 

CGCM-B23-205OS 

HADCM3-A2A-2050S 

HADCM3-B2B-2050S 

NCARPCM-A1  B-2050s 

Statistic 

1990) 

Change  (%) 

Change  (%) 

Change  (%) 

Change  (%) 

Change  (%) 

7Q1 0-Low  Flow  (m3/s) 

107 

-15% 

-13% 

-7.2% 

-2.6% 

-1.6% 

30Q  10-Low  Flow  (m3/s) 

110 

-16% 

-14% 

-7.9% 

-2.9% 

-1.9% 

NOTES: 

7Q10  -  10-year  7-day  winter  low  flow  (November/15  ~ 

April/14) 

30Q10  =  10-year  30-day  winter  low  flow  (November/15 

-  April/14) 
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Table  5.8  Effects  of  Climate  Scenarios  on  Winter  Low  Flows  in  the  Beaver  River  at  Cold  Lake 


Statistic 

Beaver  River  at  Cold  Lake  Reserve  WSC  06AD006 

Baseline  Flows  (1961- 
1990) 

CCS  RNI ES  -A1 F 1  -2020s 

CGCM-B23-2020S 

HADCM3-A2A-202QS 

HADCM3-B2B-2020S 

NCARPCM-A1B-2020S 

7Q1 0-Low  Flow  (m3/s) 

1.56 

0.874 

1.45 

1.39 

1.1' 

0 

1.24 

30Q  10-Low  Flow  (m3/s) 

1.65 

0.934 

1.53 

1.46 

1.1 

7 

1.32 

i 

Statistic 

Baseline  Flows  (1961- 

1990) 

CCSRNIES-A1  FI -2050s 

CGCM-B23-2050S 

HADCM3-A2A-2050S 

HADCM3-B2B-2050S 

NCARPCM-A1B-2050S 

7Q1 0-Low  Flow  (m3/s) 

1.56 

1.24 

0.950 

1.13 

1.0 

3 

1.06 

30Q  10-Low  Flow  (m3/s) 

1.65 

1.31 

1.00 

1.19 

1.0' 

9 

1.12 

Statistic 

Beaver  River  at  Cold  Lake  Reserve  WSC  06AD006 

Baseline  Flows  (1961- 
1990) 

CCSRNIES-A1  FI -2020s 

CGCM-B23-2G20S 

HADCM3-A2A-2020S 

HADCM3-B2B-2020S 

NCARPCM-A1  B-2020s 

Change  (%) 

Change  (%) 

Change  (%) 

Change  (%) 

Change  (%) 

7Q1 0-Low  Flow  (m3/s) 

1.56 

-44% 

|  -7.0% 

-11% 

-29° 

k> 

-20% 

30Q1 0-Low  Flow  (m3/s) 

1.65 

-44% 

-7.2% 

-12% 

-29° 

fo 

-20% 

; 

Statistic 

Baseline  Flows  (1961  - 
1990) 

CCSRNIES-A1  FI -2050s 

CGCM-B23-2050S 

HADCM3-A2A-2050S 

HADCM3-B2B-2050S 

NCARPCM-A1B-2050S 

Change  (%) 

Change  (%) 

Change  (%) 

Change  (%) 

Change  (%) 

7Q  10-Low  Flow  (m3/s) 

1.56 

-20% 

-39% 

-28% 

-34°/ 

-32% 

30Q  10-Low  Flow  (m3/s) 

1.65 

-21% 

-39% 

-28% 

-34° 

4 

-32% 

NOTES: 

7Q10  -  10-year  7-day  winter  low  flow  (November/15  ~  April/14) 
30Q10  =  10-year  30-day  winter  low  flow  (November/15  ~  April/14) 
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5.3.6  Discussion  of  Hydrologic  Effects  of  Future  Climate  Regimes 

The  results  of  the  assessment  of  the  effects  of  climate  change  on  flows  in  the  Athabasca  River  and  Beaver  River 
basins  tend  to  be  in  general  agreement  with  recent  studies  carried  out  by  independent  researchers.  Studies  on 
the  effects  of  potential  climate  change  on  flows  in  watersheds  in  Alberta  include  those  by  Kerhoven  and  Gan 
(2005)  on  the  Athabasca  River  Basin  (ARB)  and  Pietroniro  and  Toth  (2006)  on  the  South  Saskatchewan  River 
Basin  (SSRB).  Kerkhoven  and  Gan  (2006)  applied  the  MISBA  model  to  the  Athabasca  River  basin  using  the 
ERA-40  re-analysis  data  of  ECMWF  (European  Centre  for  Mid-range  Weather  Forecasts).  MISBA  is  the  ISBA 
(Interactions  Soil-Biosphere-Atmosphere  developed  by  Meteo  France)  model  modified  to  statistically  represent 
sub-grid  variability  of  soil,  vegetation  and  elevation  in  the  Athabasca  River  Basin.  MISBA  was  used  to  predict 
changes  to  the  hydrology  of  the  Athabasca  River  Basin  under  conditions  predicted  by  a  number  of  climate 
scenarios  for  the  21st  century.  Although  most  of  the  scenarios  predicted  increased  precipitation  in  the  ARB,  all 
the  scenarios  resulted  in  significantly  decreased  stream  flows  by  the  end  of  the  century  (2070-2099).  This  was 
primarily  because  of  a  predicted  decrease  in  the  size  of  the  winter  snow  pack  due  to  warmer  winters.  Warmer 
winters  result  in  less  snow  accumulation  and  increased  evaporation.  Mean  annual  flows  were  predicted  to 
decrease  by  almost  25%  by  the  last  30  years  of  the  century.  The  high  flow  season  also  became  much  shorter. 

Pietroniro  and  Toth  (2006)  present  and  discuss  the  results  of  a  study  of  the  water  availability  in  the  South 
Saskatchewan  River  Basin  (SSRB)  under  climate  change  scenarios.  The  objective  of  the  study  was  to  predict 
the  future  water  availability  in  the  SSRB  under  the  potential  impact  of  climate  change  using  hydrologic  models 
calibrated  to  SSRB  and  forced  by  downscaled  climate  scenarios  projected  by  some  selected  GCMs.  Down- 
scaled  GCMs  were  used  to  project  changes  in  local  temperature  and  precipitation  patterns.  The  climate 
information  was  then  used  to  simulate  future  river  flows  in  the  SSRB  using  the  WATFLOOD,  SACRAMENTO 
and  MISBA  hydrologic  models.  Flow  predictions  vary  by  sub-basin,  with  general  reduction  in  flows  for  the 
modelled  sub-basins  ranging  from  -13%  to  -4%. 

The  range  of  effects  predicted  by  the  HSPF  model  for  the  LARP  area  is  generally  consistent  with  the  results  of 
similar  studies  in  Alberta.  Depending  on  the  flow  statistic,  the  range  of  effects  can  span  more  than  10%,  which 
reflects  the  assumptions  inherent  in  the  formulation  of  the  climate  scenarios  as  well  as  uncertainties  in  the 
predictions.  The  IPCC  recommends  that  several  climate  scenarios  need  to  be  considered  for  assessment  of 
effects  because  a  single  climate  model  scenario  does  not  provide  a  reliable  description  of  the  climatic  evolution, 
while  ensembles  of  state-of-the-art  climate  models,  on  the  other  hand,  capture  the  main  features  of  the  past 
climatic  evolution.  Given  that  the  total  number  of  combinations  of  GCMs  and  scenarios  can  be  excessive,  a  few 
combinations  of  GCMs  and  scenarios  have  been  selected  for  this  study  so  that  the  range  of  climate  scenarios 
can  represent  dry,  median  and  wet  scenarios,  while  avoiding  excessive  runs  that  may  be  redundant  (Barrow  and 
Yu,  2005).  In  addition,  it  is  important  to  determine  the  time  horizon  of  interest  for  assessing  the  effects  of  climate 
change.  For  example,  for  water  supply  assessment,  the  planning  period  could  be  20-30  years  as  opposed  to 
100  years.  The  uncertainty  in  the  predictions  of  the  effects  of  climate  change  on  water  supply  in  the  distant 
future  may  be  larger  than  the  predicted  effects  themselves.  In  any  case,  there  is  a  greater  possibility  for 
adaptation  over  the  long-term  as  better  models  become  available  and  predictions  become  more  reliable. 

The  range  of  hydrologic  effects  predicted  by  the  FISPF  model  for  the  LARP  area  reflects  the  range  of  conditions 
that  basin  planners  should  consider.  Hence,  the  predicted  effects  using  HSPF  and  future  climate  scenarios 
provide  a  rational  basis  for  water  management  and  planning  purposes  in  the  LARP  area. 

5.4  Climate  Change  and  Climate  Variability 

While  GCMs  predict  a  progressive  increase  in  mean  temperature,  a  similar  incremental  change  in  hydro-climate 
should  not  be  expected,  however,  because  time  series  of  flows  and  climate  moisture  indices  are  dominated  by 
departures  from  mean  conditions  at  seasonal  to  multi-decadal  time  scales  (Sauchyn  et  al.  2008).  Sauchyn  et  al. 
(2008)  reports  that  quasi-periodic  cycles  dominate  the  behaviour  of  Alberta’s  hydro-climate.  Climate  variability 
(i.e.,  wet  and  dry  cycles)  is  observed  both  in  the  instrumental  (1901-2000)  and  tree-ring  (1702-2003)  record  on 
inter-annual  (2-8  years),  inter-decadal  (8-20  years)  and  multi-decadal  (20-60  years)  scales.  The  variability  has 
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5.3.6  Discussion  of  Hydrologic  Effects  of  Future  Climate  Regimes 

The  results  of  the  assessment  of  the  effects  of  climate  change  on  flows  in  the  Athabasca  River  and  Beaver  River 
basins  tend  to  be  in  general  agreement  with  recent  studies  carried  out  by  independent  researchers.  Studies  on 
the  effects  of  potential  climate  change  on  flows  in  watersheds  in  Alberta  include  those  by  Kerhoven  and  Gan 
(2005)  on  the  Athabasca  River  Basin  (ARB)  and  Pietroniro  and  Toth  (2006)  on  the  South  Saskatchewan  River 
Basin  (SSRB).  Kerkhoven  and  Gan  (2006)  applied  the  MISBA  model  to  the  Athabasca  River  basin  using  the 
ERA-40  re-analysis  data  of  ECMWF  (European  Centre  for  Mid-range  Weather  Forecasts).  MISBA  is  the  ISBA 
(Interactions  Soil-Biosphere-Atmosphere  developed  by  Meteo  France)  model  modified  to  statistically  represent 
sub-grid  variability  of  soil,  vegetation  and  elevation  in  the  Athabasca  River  Basin.  MISBA  was  used  to  predict 
changes  to  the  hydrology  of  the  Athabasca  River  Basin  under  conditions  predicted  by  a  number  of  climate 
scenarios  for  the  21st  century.  Although  most  of  the  scenarios  predicted  increased  precipitation  in  the  ARB,  all 
the  scenarios  resulted  in  significantly  decreased  stream  flows  by  the  end  of  the  century  (2070-2099).  This  was 
primarily  because  of  a  predicted  decrease  in  the  size  of  the  winter  snow  pack  due  to  warmer  winters.  Warmer 
winters  result  in  less  snow  accumulation  and  increased  evaporation.  Mean  annual  flows  were  predicted  to 
decrease  by  almost  25%  by  the  last  30  years  of  the  century.  The  high  flow  season  also  became  much  shorter. 

Pietroniro  and  Toth  (2006)  present  and  discuss  the  results  of  a  study  of  the  water  availability  in  the  South 
Saskatchewan  River  Basin  (SSRB)  under  climate  change  scenarios.  The  objective  of  the  study  was  to  predict 
the  future  water  availability  in  the  SSRB  under  the  potential  impact  of  climate  change  using  hydrologic  models 
calibrated  to  SSRB  and  forced  by  downscaled  climate  scenarios  projected  by  some  selected  GCMs.  Down- 
scaled  GCMs  were  used  to  project  changes  in  local  temperature  and  precipitation  patterns.  The  climate 
information  was  then  used  to  simulate  future  river  flows  in  the  SSRB  using  the  WATFLOOD,  SACRAMENTO 
and  MISBA  hydrologic  models.  Flow  predictions  vary  by  sub-basin,  with  general  reduction  in  flows  for  the 
modelled  sub-basins  ranging  from  -13%  to  -4%. 

The  range  of  effects  predicted  by  the  HSPF  model  for  the  LARP  area  is  generally  consistent  with  the  results  of 
similar  studies  in  Alberta.  Depending  on  the  flow  statistic,  the  range  of  effects  can  span  more  than  10%,  which 
reflects  the  assumptions  inherent  in  the  formulation  of  the  climate  scenarios  as  well  as  uncertainties  in  the 
predictions.  The  IPCC  recommends  that  several  climate  scenarios  need  to  be  considered  for  assessment  of 
effects  because  a  single  climate  model  scenario  does  not  provide  a  reliable  description  of  the  climatic  evolution, 
while  ensembles  of  state-of-the-art  climate  models,  on  the  other  hand,  capture  the  main  features  of  the  past 
climatic  evolution.  Given  that  the  total  number  of  combinations  of  GCMs  and  scenarios  can  be  excessive,  a  few 
combinations  of  GCMs  and  scenarios  have  been  selected  for  this  study  so  that  the  range  of  climate  scenarios 
can  represent  dry,  median  and  wet  scenarios,  while  avoiding  excessive  runs  that  may  be  redundant  (Barrow  and 
Yu,  2005).  In  addition,  it  is  important  to  determine  the  time  horizon  of  interest  for  assessing  the  effects  of  climate 
change.  For  example,  for  water  supply  assessment,  the  planning  period  could  be  20-30  years  as  opposed  to 
100  years.  The  uncertainty  in  the  predictions  of  the  effects  of  climate  change  on  water  supply  in  the  distant 
future  may  be  larger  than  the  predicted  effects  themselves.  In  any  case,  there  is  a  greater  possibility  for 
adaptation  over  the  long-term  as  better  models  become  available  and  predictions  become  more  reliable. 

The  range  of  hydrologic  effects  predicted  by  the  HSPF  mode!  for  the  LARP  area  reflects  the  range  of  conditions 
that  basin  planners  should  consider.  Hence,  the  predicted  effects  using  HSPF  and  future  climate  scenarios 
provide  a  rational  basis  for  water  management  and  planning  purposes  in  the  LARP  area. 

5.4  Climate  Change  and  Climate  Variability 

While  GCMs  predict  a  progressive  increase  in  mean  temperature,  a  similar  incremental  change  in  hydro-climate 
should  not  be  expected,  however,  because  time  series  of  flows  and  climate  moisture  indices  are  dominated  by 
departures  from  mean  conditions  at  seasonal  to  multi-decadal  time  scales  (Sauchyn  et  al.  2008).  Sauchyn  et  al. 
(2008)  reports  that  quasi-periodic  cycles  dominate  the  behaviour  of  Alberta’s  hydro-climate.  Climate  variability 
(i.e.,  wet  and  dry  cycles)  is  observed  both  in  the  instrumental  (1901-2000)  and  tree-ring  (1702-2003)  record  on 
inter-annual  (2-8  years),  inter-decadal  (8-20  years)  and  multi-decadal  (20-60  years)  scales.  The  variability  has 
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been  linked  to  sea  surface  temperature  (SST)  anomalies  in  the  tropical  and  extra-tropical  regions  of  the  Pacific 
Ocean,  ENSO  and  PDO,  respectively. 

The  Prairies  have  Canada’s  most  variable  hydro-climate,  and  the  wet  and  dry  cycles  have  significant 
implications  for  prairie  communities  and  economies  and  the  management  of  land  and  water  resources. 
Therefore,  understanding  the  nature  and  causes  of  this  inter-annual  to  multi-decadal  variability  in  hydro-climate 
is  important  for  an  integrated  assessment  of  the  effects  of  climate  variability  on  water  yield  versus  those  due  to 
climate  change.  Sauchyn  et  al.  (2008)  recommend  that  an  additional  criterion  for  evaluating  the  applicability  of 
GCMs  should  be  their  ability  to  simulate  inter-annual  to  decadal  climate  variability,  particularly  for  a  climate  like 
Alberta’s. 

The  calibrated  HSPF  model  calibrated  for  the  Athabasca  River  Basin  was  used  on  the  Firebag  River  sub-basin 
to  generate  a  baseline  (1961-1990)  as  well  as  a  55-year  (1954-2007)  series  of  flows  using  available  data  at  the 
Fort  McMurray  Airport  climate  station.  Figure  5.5  shows  that  the  observed  and  simulated  flow  statistics  for  the 
period  of  record  are  comparable.  Appendix  K  provides  detailed  information  on  the  analysis  carried  out  on  the 
Firebag  River  sub-basin. 

It  is  assumed  that  the  simulated  55-year  flow  data  series  may  capture  the  variability  due  to  past  ENSO  events 
and  possibly  some  past  PDO  events  to  a  greater  extent  than  the  standard  30-year  (1961-1990)  series.  Table 
5.9  shows  the  forecasted  changes  in  flow  statistics  for  the  1954-2007  period  due  to  the  2020s  and  2050s  climate 
change  scenarios. 

The  effect  of  future  climate  scenarios  on  flow  variability  was  assessed  by  comparing  the  standard  deviation  and 
coefficient  of  variation  for  the  two  (1961-1990  and  1954-2007)  annual  flow  series  under  the  baseline,  2020s  and 
2050s  climate  scenarios.  Figure  5.6  shows  the  simulated  annual  mean  flows  for  the  1954-2007  period  when 
observed  climate  data  are  available  and  the  flow  series  using  the  2020s  and  2050s  climate  change  scenarios.  It 
is  apparent  from  Figure  5.6  that  the  period  from  1954  to  1974  was  relatively  wet  compared  to  the  period  from 
1975  to  2007,  with  the  period  from  1954  to  2007  likely  reflecting  a  low  frequency  wet-dry  cycle.  The  period  from 
1961  to  1990  (standard  climate  baseline  period)  would  also  encompass  portions  of  the  wet  and  dry  cycles. 
Table  5.10  compares  the  statistics  of  the  1961-1990  and  1954-2007  series  using  the  2020s  and  2050s  climate 
change  scenarios.  The  coefficients  of  variation  shown  in  Table  5.10  suggest  that  there  is  no  difference 
(difference  not  statistically  significant)  between  the  baseline  1961-1990  and  1954-2007  series,  as  well  between 
the  same  two  series  under  the  five  climate  scenarios  for  the  2020s  and  2050s.  One  reason  for  the  lack  of 
difference  between  the  two  series  (30-year  and  54-year)  is  likely  because  the  54-year  is  still  too  short  and  does 
not  include  the  extreme  wet  and  dry  cycles  of  the  past  as  suggested  by  analysis  of  tree-ring  records.  It  is 
important  to  consider  the  probability  and  potential  for  the  occurrence  of  very  extreme  events,  either  excessive 
moisture  or  drought,  that  are  not  reflected  in  the  30-yearl 961 -1990  baseline  for  “current  climate”.  However, 
since  the  interaction  of  the  various  natural  climate  cycles  is  not  entirely  understood,  and  the  complicating  effects 
of  a  warmer  climate  add  another  layer  of  complexity,  this  is  a  relatively  new  area  of  research  in  climate  science. 
Tree  ring  records,  for  example,  indicate  that  longer,  more  frequent  and  more  severe  drought  periods  have 
occurred  in  the  past  on  the  prairies,  but  linking  climate  change  feedbacks  to  these  observations  to  express  the 
probability  and  severity  of  these  extreme  events  in  the  future  is  very  challenging. 

The  assessment  of  climate  variability  in  the  LARP  area  will  always  be  constrained  by  the  lack  of  relatively  long 
(>  100  years)  instrumental  data  records.  Analysis  of  tree-rings  can  provide  some  information  on  past  climate 
and  streamflow  variability,  however,  this  information  cannot  be  replicated  with  or  incorporated  into  continuous 
simulation  hydrologic  models  such  as  HSPF.  One  alternative  approach  is  to  generate  synthetic  series  of 
temperature,  precipitation  and  flows  using  the  statistical  characteristics  of  past  and  present  tree-ring  records  in 
conjunction  with  the  relationship  between  present  tree-ring  records  and  observed  climate  and  flow  data  series. 
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5.5  Summary 

To  address  the  potential  effects  of  climate  change  on  water  yield  in  the  Lower  Athabasca  Regional  Plan  Area 
(LARP),  which  includes  the  downstream  portion  of  the  Athabasca  River  Basin  and  the  Beaver  River  Basin,  the 
HSPF  model,  calibrated  and  validated  for  these  two  basins,  has  been  used  to  simulate  the  hydrologic  effects  of 
forecasted  future  climate  scenarios.  The  1961  to  1990  period  was  used  as  the  climatological  baseline  period  for 
the  effects  assessments. 

AENV  provided  Golder  with  the  Alberta  Climate  Model  data  and  average  monthly  changes  in  temperature  and 
precipitation  predicted  by  five  GCMs  for  two  future  periods:  2010  to  2039  (referred  to  as  2020s)  and  2040  to 
2069  (referred  to  as  2050s).  The  differences  between  the  baseline  averages  and  the  2020s  averages,  and  the 
differences  between  the  baseline  averages  and  the  2050s  averages  were  estimated.  The  average  change  in 
mean  annual  precipitation  varies  from  -1%  to  +7%  for  the  2020s  scenarios  and  from  +2%  to  +13%  for  the  2050s 
scenarios.  The  range  of  the  change  in  precipitation  is  much  wider  on  a  monthly  basis:  -10%  to  +23%  for  the 
2020s  scenarios  and  -15%  to  +46%  for  the  2050s  scenarios.  The  average  increase  in  mean  annual  temperature 
varies  from  0.63°C  to  1 ,52°C  for  the  2020s  scenarios  and  from  1 ,77°C  to  4.35°C  for  the  2050s  scenarios. 

The  HSPF  model  calibrated  for  the  Athabasca  River  Basin  and  Beaver  River  Basin  was  run  with  the  baseline 
climate  data  and  the  adjusted  future  climate  data.  The  effects  on  flows  were  assessed  using  the  flow  statistics  at 
selected  locations  on  the  main  stem  of  the  Athabasca  River  and  on  the  Beaver  River  at  the  Cold  Lake  Reserve 
for  the  2020s  and  2050s  climate  scenarios. 

General  conclusions  are  as  follows: 

■  All  the  five  climate  scenarios  for  the  2020s  tend  to  result  in  lower  mean  annual  flows,  with  the  decrease  in 
mean  annual  flow  becoming  more  severe  downstream  along  the  Athabasca  River. 

■  Changes  in  mean  monthly  flows  for  the  2020s  scenarios  tend  to  be  positive  (increase  relative  to  baseline 
values)  during  the  spring  months  and  significantly  negative  (decrease)  during  the  summer  months.  The 
seasonal  differences  in  flows  are  consistent  with  the  relative  effects  of  increased  precipitation  (more  spring 
runoff)  and  increased  temperatures  during  the  summer  months  (higher  summer  evapotranspiration). 

■  The  2050s  climate  scenarios  result  in  changes  in  mean  annual  flow  that  are  generally  more  severe 
compared  to  the  2020s  scenarios,  which  is  consistent  with  the  forecasted  changes  in  temperature  and 
precipitation. 

■  The  range  of  effects  of  future  climate  regimes  on  flows  in  the  Lower  Athabasca  Regional  Plan  Area  (LARP) 
is  summarized  below.  Three  locations  are  used  for  the  summary:  (1)  Athabasca  River  at  Athabasca 
(07BE001)  to  represent  inflow  to  the  LARP  area,  (2)  Athabasca  River  below  McMurray  (07DA001)  to 
represent  flows  in  the  oil  sands  region  where  the  amount  of  water  withdrawal  is  a  key  consideration  in  the 
planning  process,  and  (3)  Beaver  River  at  Cold  Lake  Reserve  sub-basin  (06AD006)  to  represent  the 
southern  portion  of  the  LARP  area. 

The  summary  results  indicate  that  water  yield  from  the  two  basins  will  generally  decrease,  with  the  effects 
(in  percentage  terms)  more  significant  during  the  month  of  August.  The  effects  on  low  February  flows, 
although  reduced,  tend  to  be  less  (in  percentage  terms)  compared  to  changes  in  August  flows.  Flood  flows 
tend  to  be  higher  under  the  future  climate  scenarios. 

The  results  of  the  assessment  of  the  effects  of  climate  change  on  flows  in  the  Athabasca  River  and  Beaver  River 
basins  tend  to  be  in  general  agreement  with  recent  studies  carried  out  by  independent  researchers.  Kerkhoven 
and  Gan  (2006)  predicted  that  mean  annual  flows  in  the  ARB  would  decrease  by  almost  25%  by  the  last 
30  years  of  the  century.  Pietroniro  and  Toth  (2006)  predict  general  reduction  in  flows  in  the  SSRB  ranging  from  - 
13%  to -4%. 
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The  range  of  effects  predicted  by  the  HSPF  model  for  the  LARP  area  is  generally  consistent  with  the  results  of 
similar  studies  in  Alberta.  Depending  on  the  flow  statistic,  the  range  of  effects  can  span  more  than  10%,  which 
reflects  the  assumptions  inherent  in  the  formulation  of  the  climate  scenarios  as  well  as  uncertainties  in  the 
predictions.  The  range  of  hydrologic  effects  predicted  by  the  HSPF  model  for  the  LARP  area  reflects  the  range 
of  conditions  that  basin  planners  should  consider.  Hence,  the  predicted  effects  using  HSPF  and  future  climate 
scenarios  provide  a  rational  basis  for  water  management  and  planning  purposes  in  the  LARP  area. 

The  scope  of  work  for  this  study  included  an  assessment  of  flow  variability  relative  to  effects  due  to  forecasted 
climate  scenarios.  The  calibrated  HSPF  model  calibrated  for  the  Athabasca  River  Basin  was  used  on  the 
Firebag  River  sub-basin  to  generate  a  baseline  (1961-1990)  as  well  as  a  54-year  (1954-2007)  series  of  flows 
using  available  data  at  the  Fort  McMurray  Airport  climate  station.  It  was  assumed  that  the  simulated  55-year 
flow  data  series  may  capture  the  variability  due  to  past  ENSO  events  and  possibly  some  past  PDO  events  to  a 
greater  extent  than  the  standard  30-year  (1961-1990)  series.  The  effect  of  future  climate  scenarios  on  flow 
variability  was  assessed  by  comparing  the  standard  deviation  and  coefficient  of  variation  for  the  two  (1961-1990 
and  1954-2007)  annual  flow  series  under  the  baseline,  2020s  and  2050s  climate  scenarios.  The  calculated 
coefficients  of  variation  suggest  that  there  is  no  difference  (difference  not  statistically  significant)  between  the 
baseline  1961-1990  and  1954-2007  series,  as  well  between  the  same  two  series  under  the  five  climate 
scenarios  for  the  2020s  and  2050s.  One  reason  for  the  lack  of  difference  between  the  two  series  (30-year  and 
54-year)  is  likely  because  the  54-year  is  still  too  short  and  does  not  include  the  extreme  wet  and  dry  cycles  of  the 
past  as  suggested  by  analysis  of  tree-ring  records. 
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6.0  CLOSURE 

This  report  presents  the  outcome  of  the  assessment  of  hydrologic  models  for  the  Lower  Athabasca  Regional 
Plan  Area  and  the  results  of  the  HSPF  model  calibration  and  validation,  and  simulations  with  baseline  (1961- 
1990)  and  climate  change  (2020s  and  2050s)  scenarios.  Please  direct  any  questions  or  clarifications  regarding 
the  contents  of  this  draft  report  to  Anil  Beersing  at  (403)  260  2292. 
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7.0  THIRD  PARTY  DISCLAIMER 

This  draft  report  has  been  prepared  by  Golder  Associates  Ltd.  (Golder)  for  the  benefit  of  the  client  to  whom  it  is 
addressed.  The  information  and  data  contained  herein  represent  Golder's  best  professional  judgment  in  light  of 
the  knowledge  and  information  available  to  Golder  at  the  time  of  preparation.  Except  as  required  by  law,  this 
draft  report  and  the  information  and  data  contained  herein  are  to  be  treated  as  confidential  and  may  be  used  and 
relied  upon  only  by  the  client,  its  officers  and  employees.  Golder  denies  any  liability  whatsoever  to  other  parties 
who  may  obtain  access  to  this  draft  report  for  any  injury,  loss  or  damage  suffered  by  such  parties  arising  from 
their  use  of,  or  reliance  upon,  this  draft  report  or  any  of  its  contents  without  the  express  written  consent  of  Golder 
and  the  client. 
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Table  A.3  Assessment  of  MISBA  for  Hydro-Climate  Modelling  of  Lower  Athabasca  Region 


t  of  MIKE-SHE  for  Hydro-Climate  Modelling  of  Lower  Athabasca  Region 


Table  A.5  Assessment  of  HSPF  for  Hydro-Climate  Modelling  of  Lower  Athabasca  Region 


Table  A.7  Assessment  of  HYDROTEL  for  Hydro-Clin 


Table  A.8  Assessment  of  SSARR  for  Hydro-Climate  Modelling  of  Lower  Athabasca  Region 


Sub-Criterion 

— 

event  to  Objectives  of  Study 

How  important  is  this 

SSARR 

1 - - - - 1 

Ability  to  Incorporate  Key  Features  of  Watershed  being  Modelled 

_ 

i  i 

1 

Physical  Characteristics  such  as  Area,  Elevation,  Slope 

The  model  considers  drainage  area  but  not  slope.  The  catchment  can  be  divided  into  elevation  band,  and  elevation  is  used  to  represent  bands  in  a  snow  melt  model. 

: 

Soil/Surficial  Geology 

|Not  considered 

l 

3 

Veqetation  Type/Density 

|Not  directly,  but  has  a  simple  function  that  simulates  the  interception  process  of  vegetation.  A  maximum  interception  quantity  based  on  the  type  and  density  of  veqetation  can  be  specified 

1  i 

4 

Aspect 

Not  considered 

5 

Anthropogenic  Changes  in  Land  Cover 

Not  considered 

6 

Variable  Runoff  Contributing  Areas  (e.g.  Non-contributinq  Areas  under 

Not  considered 

Ability  to  Simulate  Dominant/Key  Hydrologic  Processes  in  Watershed  being  Modelled 

1 

i 

7 

Direct  Runoff  from  Rain 

Runoff  from  rainfall  is  calulated  from  empirically  derived  relationships  of  soil  moisture  index  versus  runoff  percent.  Input  intensity  can  be  included  as  a  third  variable. 

8 

Direct  Runoff  from  Snowmelt 

Snowmelt  is  calculated  by  the  temperature  index  or  energy  budget  approach.  The  model  also  distinguishes  melt  during  clear  weather  and  rain  and  account  snowpack  conditioninq  which  accounts  for  the  "cold  content"  and  liaui 

d  water  deficiency  of  the  snowpack. 

9 

Interflow 

The  runoff  has  surface  and  subsurface  or  interflow  components. 

10 

nfiltration/Groundwater 

Not  explicitly  simulated;  implicitly  incorporated  into  baselow. 

11 

nterception 

Maximum  interception  quantity  based  on  the  type  and  density  of  vegetation  is  specified  by  the  user.  A  simple  function  in  the  model  simulates  the  interception  process 

12 

3asef!ow 

The  base  flow  components  are  computed  as  the  product  of  Base  Flow  Percent  and  Runoff  Rate.  Runoff  that  contributes  to  base  flow  is  estimated  as  a  function  of  the  base  flow  infiltration  index,  an  index  of  depression  sforaqe. 

Miich  holds  runoff  available  for  deep  percolation. 

13 

Winter  flows 

The  base  flow  component  has  a  longer  term  and  shorter  term  components.  This  separation  of  components  gives  the  watershed  model  the  capability  of  reproducing  the  late  summer,  fall  and  winter  qroundwater  recession. 

14 

Lake  Evaporation 

'Jot  considered 

15 

Evapotranspiration 

Potential  Evapotranspiration  (PET)  is  computed  using  Thornthwaite’s  formula  converted  to  produce  daily  values.  A  temperature  vs.  Adjustments  for  season  and  latitude,  elevation  and  rainy  days  are  made. 

16 

Snow  sublimation 

'Jot  considered 

17 

Snow  transport/re-distribution 

Not  directly,  but  adjustment  can  be  made  for  gauge-catch  deficiency  before  calculating  the  basin  weighted  average  precipitation,  through  the  user-specified  characteristics  of  the  precipitation  gauge. 

Ability  to  Simulate  Processes  on  Event  to  Multi-year  Basis 

18|Event-simulation  only 

Yes 

19|Continuous  Simulation 

Yes 

Abitity  to  Simulate  Key  Hydraulic  Features 

20 

Main  Channel  Routing 

Yes.  The  program  allows  flexibility  to  determine  routing  coefficients  which  simulate  downstream  peaks  and  timing  response.  Channel  routing  can  be  accomplished  with  either  a  routing  equation  for  incremental  routing,  or  a  table  which  specifies  time  of  storage-discharge  relationships. 

21 

Routing  through  a  Network  of  Channels 

Yes.  The  program  allows  flexibility  to  determine  routing  coefficients  which  simulate  downstream  peaks  and  timing  response.  Channel  routing  can  be  accomplished  with  either  a  routing  equation  for  incremental  routing,  or  a  table  which  specifies  time  of  storage-discharge  relationships. 

22 

Storage  in  Lakes/Wetlands 

Yes.  Routing  of  natural  lake  is  based  on  free-flow  conditions,  i.e.,  elevation-outflow  relationships  are  fixed,  and  outflow  is  determined  by  hydraulic  head.  Routing  is  accomplished  by  an  iterative  solution  of  the  continuity  storage  equation. 

23 

=loodplain  Routing  and  Storage 

Yes,  through  a  "cascade  of  reservoirs"  technique,  wherein  the  lag  and  attenuation  of  the  flood  wave  are  simulated  through  successive  increments  of  lake-type  storage. 

24 

-low  under  Ice  Cover 

Yes.  The  snow  pack  may  be  melted  by  snow-ground  interface  and  groundmelt  is  accounted  for  by  a  month  vs.  groundmelt  relation,  and  the  resulting  water  enters  the  runoff  simulation  without  being  affected  by  the  snow  conditioninq. 

25 

ce  jams  and/or  break-up 

Not  considered 

Abddy  to  Provide  Outputs  a.  Spatial  Scale  o,  In.eree.  within  Watershed  being  Modelled 

! 

26 

Watershed 

Yes 

27 

Sub-basin 

Yes 

28 

Local  Drainage/Landscape  Feature 

Yes 

29 

Hillslope  catchment 

No 

30 

Lakes/Wetlands 

No 

31 

Perched  Basins 

No 

32 

Beaver  Dams 

No 

Ab®ty  to  Provide  Outputs  at  Time  Scale  of  Interest 

1 

33|Annual 

Yes 

34  Monthly 

Yes 

35Daily 

Yes 

36 1  Hourly 

No 

AbiSty  of  Model  to  Use  Format  of  Downscaled  Cimate  Data 

i 

37|With  no  or  minimal  pre-processing 

No 

|With  considerable  pre-processing 

Yes 

Climate 

391 

Point 

1-hour  to  1-day  average 

Grid  cell 

No 

Sub-basin 

1-hour  to  1-day  average 

- -  .  , 

Point 

1-hour  to  1-day  average 

Grid  cell 

No 

Sub-basin 

1-hour  to  1-day  average 

Isohyets 

No 

Thiessen  Polygons 

- - L . . . — - - ! 

Snow  pillow 

Average  snow  depth 

No 

Snow  depth  by 

No 

Snow  density 

iNo  ’  - - 

Snow  water  equivalent 

42 

Hourly 

Yes 

Daily 

Yes 

43 

Hourly 

Daily 

Yes 

44 

Relative  Humidity 

Daily 

Yes.  It  is  calculated  from  Temperature  and  Dew  Point  Temperature 

45 

Cloudiness 

Daily 

Yes,  through  Insolation 

46 

Evaporation 

Daily 

47 

Evapotranspiration 

Daily 

Yes 

Gnd  cells 

L 

No 

4S 

Sub-basins 

Yes 

5C 

)  Elevation  contours 

Yes 

51 

1  DEM 

No 

- - - - | — — — | - 

5: 

l  GIS-based  Hydrography 

V  '  ^  ‘  " - - - 1 - 1 - - 1 

J  Conceptual/ Schematized  Network  | 

Yes 

Criterion 

Sub-Criterion 

54 

Storage  Elements 

Yes 

I - - - - - - - - - - - - - 

55 

56 

Surficial 

No 

Soil 

No  - - - - „■  ^  ,  -  , - 

IVaqatafion  Coverage 

57 

No 

58 

Forest 

Not  directly  but  interception  is  quantified  by  the  user  based  on  the  type  and  density  of  vegetation. 

59 

Grass 

No 

60 

No 

Criteria  and  Sub-Criteria  for 

Degree  of  Sophistication  of  Processes  Simulated  In  Model  - 

61 

Unit  Hydrograph 

No 

62 

Conceptual  storages 

Yes 

63 

Hydrologic  Response  Units 

No 

64 

Water  Balance 

No 

1 - - - = - - - - - - - - 

65 

Degree  day 

Yes  1.11 

66 

Energy  Equations 

Yes 

67 

Energy  Equations 

No 

68 

Morton's  Approach 

No.  Uses  Thomthwaite's  formula 

69 

Pan  Evaporation 

Mo 

I 

70 

Hydrograph  Translation 

Yes 

71 

Muskingum 

Yes 

72 

Physics-based:  Spatial  and  Temporal  Distribution 

No 

73 

Average  based  on  soil  type 

No 

Glacier-melt 

! 

74 

Based  on  temperature 

No 

75 

Steady  rate 

No 

Criteria  Defining  Ease  of  Us 

by  Modeler 

1 - 

76 

Commercial-Wide  Use 

L . 

77 

Commercial-Limited  Use 

78 

3ublic  Domain-Wide  Use 

Yes 

79 

3ublic  Domain-Limited  Use 

80 

Research-Real  World  Applications 

81 

Research-Limited  Application 

82 

Research-Under  Development 

Traning/Support 

83 

Training  provided  by  developers 

84 

Training  provided  by  consultants 

85 

Frequent  training  workshops 

86 

Help  Files  and  Example  Applications 

87 

Win  95 

Yes 

88 

Win  98 

Yes 

89 

Win  2000 

Yes 

90 

UNIX 

91 

DOS 
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